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An experimental and analytical investigation was conducted to 
determine the free surface shapes of circular liquid jets impinging 
normal to sharp-edged disks under both normal and zero gravity con- 
ditions. An order of magnitude analysis was conducted indicating re- 
gions where viscous forces were not significant when computing free 
surface shapes. The demarcation between the viscous and Inviscid re- 
gion was found to depend upon the flow Reynolds number and the ratio 
between the jet and disk radius. 

Experiments conducted under zero gravity conditions yielded 
three distinct flow patterns. These flow patterns were defined as 
surface tension flow, transition flow, and inertia flow. The flow 
regions were classified in terms of the relative effects of surface 
tension and inertial forces. The transition between regions was cor- 
related with the system Weber number and the ratio of the jet to the 
disk radius. The normal gravity plume shapes were observed to jump 
from one apparently stable flow pattern to another until steady-state 
was reached. 

A zero gravity inviscid analysis was performed in which the 
governing equations and boundary conditions in the physical plane 
were transformed into an inverse plane. In the Inverse plane, the 
stream function and velocity potential became the coordinates thus 
removing the prime difficulty In free surface problems, that of hav- 
ing to guess at the true position of the free surface. The governing 
equations were nonlinear In the Inverse plane thus requiring a nu- 
merical solution in which sets of nonlinear algebraic equations were 
solved simultaneously. Comparisons between experiment and numerical 
computations were made for the infinite and finite plate cases with 
the result that good agreement for the free surface shapes wore ob- 
tained. 
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I. INTRODUCTION 


A krawXedge of the dynamics of free liquid jets is required for 
the solution of a variety of problems associated with fluid flow 
within propellant tanks under low gravitational conditions. In par- 
ticular, an understanding of the liquid Jet - Impact process, as 
occurs when liquid Impinges upon baffles or tank walls during an in- 
flow or reorientation maneuver, will be required In order to predict 
llquld-propellant location, heat transfer rates and pressure distribu- 
tions. The area of liquid jet Impingement also has direct applicabil- 
ity to the spacecraft fire safety problem, In which water jets are 
employed as extlngulshant agents under low gravity conditions. In 
order to be 'able to predict required delivery flow rates, the accurate 
prediction of flow surface coverage as a function of jet momentum is 
needed . 

There generally appears to be three chief obstacles which have 
in the past prohibited the attainment of solutions to steady-state 
liquid jet-solid interaction problems. The major obstacle is the 
presence of the free surface. In order to apply numerical techniques 
to the solution of free-jst problems it is necessary to define the area 
over which the computations are made by boundaries defined by the free 
liquid surface. Unfortunately, the location of the free-surface is 
one of the items sought from the solution so that various techniques 
must be devised to circumvent this situation. Furthermore, analytical 
techniques are restricted solely to two-dimensional problems, whether 
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a free surface exists or not. The second obstacle is gravity. Liquid 
Jets In air (free jets), unlike llquld-into-llquld jets and gas-into- 
gas Jets (submerged Jets), are affected significantly by gravitational 
forces. The free~surface shape and velocity profiles are dependent on 
both the magnitude and the orientation of gravity. The addition of 
gravity necessarily complicates a model either through the governing 
equation or through the boundary conditions. Neglecting gravity In 
the model makes questionable the comparison of the theory with normal- 
gravity experimental data. The final obstacle Is surface tension, an 
effect which has generally been neglected in almost all studies on 
free Jets. The addition of surface tension Into a model leads to non- 
linear free surface boundary conditions. 

The purpose of this report Is to present the results of an exper- 
imental and analytical study conducted at the NASA Lewis Research Cen- 
ter concerning zero gravity Isothermal liquid Jet Impingement. An 
axis 'mmetrlc liquid Jet was impinged normally onto a sharp-edged disk 
under conditions In which both inertial and surface tension forces 
are of Importance. The experimental free surface shapes were corre- 
lated with known system parameters. An analytical model was formu- 
lated and the free surface shapes and streamlines were calculated for 
a number of discrete cases. 
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II. LITERATURE SURVEY 
A. Experimental Studies 

Very few experimental studies have been conducted to examine free 
jets impinging on solid surfaces. No work has been conducted where 
the major concern was either the shape of the free surface or the meas- 
urement of velocity profiles within the jet. Also, only one experi- 
ment has been conducted using a two-dimensional jet. A two-dimensional 
jet is one in which the flow emanates from a rectangular slot in which 
the width of the jet is very large relative to the thickness of the jet. 
Schach (38) measured the pressure distribution and analytically calcu- 
lated the free surface shape and velocity distributions for jets im- 
pinging onto flat panels at various impingement inclinations relative 
to the direction of flow. The jet employed had dimensions of 21 by 
115 mm. According to Schach, the jet diverged spatially after imping- 
ing upon the panel and thus it can only be considered as truly two- 
dimensional close tu the centerline. An excellent account of an 
elaborate experimental apparatus for obtaining a quiescent circular 
water jet in normal gravity is given by Donnelly, et al. (11). 

Their major concern was jet stability under imposed audio frequency 
disturbances and, therefore, the impingement phenomenon was not di- 
rectly observed. Rupe (37) and Stephens (41) experimentally meas- 
ured the pressure distribution caused by circular jets striking solid 
surfaces in normal gravity. However, neither Rupc or Stephens meas- 
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ured the free-surface shape or discussed any instabilities which oc- 
curred. 

In nearly all flows where a circular liquid Jet strikes a large 
flat surface, typically what happens in normal gravity is that the 
liquid Jet impinges on the surface and moves radially outward from 
the stagnation point until a certain radial distance is reached 
whereupon an instability known as a circular hydraulic jump occurs. 
The Jump is characterized by an abrupt Increase in the liquid depth 
and turbulent fluid motion. Koloaeua, et al, (22) were concerned 
solely with predicting the behavior of the circular hydraulic jump. 

A water Jet impinging on a flat plate of epoxy material was employed 
in these experiments. The circular hydraulic Jump was the subject 
for a very complete study conducted by Nirapathdongporn (33) , whose 
report contains an excellent description of various devices for meas- 
uring Jet shapes and Jet diameters. 

All of the above mentioned studies deal with normal-gravity, 
liquid Jet-solid impingement. There have been no experimental studies 
on the impingement of liquid Jets under zero-gravity conditions. 

B. Analytical Studies 

!• Steady two-dimensional potential flow . - A number of papers 
and books have discussed steady-state two-dimensional free jets im- 
pinging on a variety of surfaces using analytical techniques. The 
majority of these studies wore concerned with Irrotational, incom- 
pressible, Invlscid flow, in which the effects of gravity and surface 
tension were neglected. One of the major attractions of this type of 
problem la that It can be handled using complex potential theory and, 
therefore, can be treated analytically. 
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A two-dimensional jet striking an infinitely flat surface at 
various angles was examined by Batchelor (4), who solved for the limit- 
ing stream thickness as a function of flow Impingement angle and jet 
diameter. However, no attempt was made to predict free-streamline 
shape. Schach (38) treated the impingement as a function of angle 
using Prandtl's hodograph method, and obtained the equations for the 
free surface shapes, flow distribution, and pressure distribution for 
the case of Impingement on an infinitely wide plate. Kochin, et 
al. (21) also examined the. impingement of a two-dimensional jet 
obliquely to an infinite flat plate, and discussed the case of im- 
pingement on a plate of finite length. The equation of the free- 
surface for the case of a two-dimensional jet striking « flat surface 
at right angles is presented by Mllne-Thomson (30) who also solved 
for the velocity components within the jet. An excellent discussion 
of the techniques for handling two-dimensional free jet problems is 
presented by Gurevich (14) . Some of the two-dimensional flows ex- 
amined by Gurevich Include flow around a finite wedge, perpendicular 
to a finite plate, obliquely to an infinite fi'»t plat-, and flows 
where a variety of solid objects are positioned adjacent to one wall 
or between two walls. Chang, et al. (9) analyzed the two-dimensional 
flow of a jet interacting with a number of flat segments at angles to 
one another. The results include flow turning angles but not free 
surface shapes or velocity profiles. The Irrotational flow pattern 
of a free jet discharging from a slot and flowing past a wedge was 
analyzed by Arbhabhlrama (3). All of the above texts and articles 
were concerned with analytical techniques for obtaining solutions. 

The area of steady two-dimensional potential flew represents tiie most 
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complete area of research in the field of jet Impingement, 

Steady axlaymmetrlc vlscoua How. - Watscn ( 43 ) has analyti- 
cally investigated free Jet-lmpingemenL for the case of large Reynolds 
numbers where the viscous forces are confined to a thin boundary layer 
adjacent to the plate. A slmiUrlty solution was obtained for both 
the two-dimensional and axisymmetric velocity profiles and free sur- 
face shapes for the case of normal impingement. As mentioned by 
Watson, the similarity solution can only be expected to be valid when 
the radial distance is sufficiently large for the Incident jet to have 
lost its influence. The effects of gravity and surface tension were 
neglected in the analysis. Watson solved for the radial position of 
the circular hydraulic Jump. 

C. Numerical Studies 

Steady two-dimensional potential flow . - When the shape of 
the solid upon which the jet Impinges becomes complex, numerical 
techniques for the solution of free jet problems have to be applied. 
Jeppson (19) presents an excellent article in this regard. Jeppson 
employed the stream function and the velocity potential as the inde- 
pendent variables and the coordinates as the dependent variables. A 
similar inversion approach has been previously used to solve a variety 
of fluid dynamics problems as shown in references 5, 20, 31, 42, 
and 44 and is mainly attributable to Thom and Apelt (42) . Using this 
technique, Jeppson was able to circumvent the problem of working In 
the physical plane and having to guess ut the true position of the 
free surface. The latter iterative approach was used in references 
1» 8, 13, 27, 32, 36, and 40 with limited success. Jeppson solved 
the problem the two-dimensional flow over a wedge, and, as such, la 
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the only one to have attempted numerical solutions of this problem. 
Lastly, Chan (7) applied the finite element method to a number of 
free-Burface flow problems, including the flow from a circular ori- 
fice. 

2. Steady axisvmmetric potential flow . - The solution of axi- 
symmetric flow problems cannot utilize the powerful tool of complex 
analysis. For this reason, only numerical solutions can be attempted 
for problems of this nature. 

LeClerc (25) studied the in^ingement of an axially symmetric 
liquid jet perpendicular to a flat surface. The shape of the free 
surface was found using an electrical analogy. This method thus 
fixed the position of the free surface and enabled the author to ap- 
ply standard finite-differencing methods and employ Southwell’s re- 
laxation technique to solve Laplace’s equation. Jeppson (19) applied 
his inversion technique to find the flow pattern and free-surface 
shape for the case of axisymmetric flow past a variety of bodies of 
revolution, including cones. Jeppson also applied his technique to 
the solution of a jet of Invlscid, incompressible fluid issuing from 
a nozzle into the free atmosphere. He indicates how his method may 
be extended to a variety of other problems. Schach (39) used a semi- 
analytical technique based on Trefftz’s approximate method to find 
the shape of an axisymroetrlc free jet impinging normally on a plate. 
Also presented in Schach's article was the pressure distribution on 
the plate which was calculated from the velocity distribution using 
Bernoulli's equation. Young, et al. (45) and Brvnauer (6) determined 
the flow pattern past two disks immersed in axlsymmetrlc flow. Both 
Young and Brunauer solved Laplace’s equation in the physical plane. 
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No analyses have been conducted for the case In which an Invlacld 
free Jet impinges upon a plate of finite thickness. 

The articles mentioned above encompass all the known solutions 
with regard to axlsymnietrlc jet impingement. References 7, 8, 20, 
27, and 40 deal specifically with numerical methods applied to free 
surface problems In which no impingement occurs. Jeppson (20) em- 
ployed an inverse formulation while the others worked In the physi- 
cal plane. 

Unsteady two-di mensional and axlsvmmetrlc potential flow . - 
Huang (17, 18) has investigated unsteady flows and considered the 
impact phenomena for both two-dimensional and axisymmetric jets. 

The major interest in these articles was In obtaining the Initial 
pressure distribution due to liquid impact. 

Steady potential flow including gravitational effects . - 
The addition of gravity in analyses for potential flows causes no 
serious formulation problem for either the two-dimenslonai or axi- 
symmetric case. The reason for this is because its effect enters 
only through the free-surface boundary conditions and not the gov- 
erning equations. Jeppson (19) included gravity in his analysis of 
the impingement on a two-dimensional wedge. Moayerl, et al. (31) 
Southwell, et al. (40) and Chan (8), all considered the effect of 
gravity in dealing with steady, potential, free-surface problems in 
which no impingement occurs. 

5. ^tea^ potential f jow including surface tension . - Zhukovskll 
(46) has indicated how to Include the effects of surface tension. He 
examined a two-dimensional problem using complex analysis, but his 
method 1s not extendable to either axisymmetric or three-dimensional 
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flows. 

6, St eady thrcc-dJ mcnBlonoT potential flow. - Until v^ry recent- 
ly, very little had been accomplished in the area of three-dimensional 
potential flow with a free surface, much leas including impingement. 
Davis and Jeppaon (10) developed a computer program to solve free- 
surface problems of this type using the inverse method. Michelson 
(28, 29) also examined jets under these conditions. He treated the 
case of an axisymmetric jet Impinging obliquely on a flat surface, 
and analytically showed the occurrence of wedge-shaped dry stones when 
the impingement angle was less than a critical value. Free-surface 
shapes are not obtainable using Michelson 's method. 
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III. ORDER OF MAGNITUDE AHAI.YSIS 

A. Formulation 

The problem under conalderetlon Is the viscous flow of a circular 
liquid jet as It impinge, normall, to an Infinite flat plate, as shown 
in Plgure 1. The objective la to determine the free surface shape of 
the impinging liquid and the velocity profiles within the Jet. In gen- 
eral, flows of the type described will depend on viscous, surface ten- 
sion, inertial and body forces. Physical Intuition tells us that If 
the velocity Is large and the diameter of the plate Is sufficiently 
small, there will be regions wherein viscous forces are not of prime 
importance In determining the resulting flow behavior, particularly 
the free surface shape. The dscous forces. In this case, will be con- 
fined to a thin boundary layer on the plate which originates fro. the 
stagnation point. The location of the stagnation point Is shown In 
Figure 1. The Jet or nossl. radius Is R„ and the distance between 
fhe plate and nossle Is given as H. A cylindrical coordinate system 
(r.s) emanating from the stagnation point Is chosen. An order of mag- 
nitude analysis will permit the governing equations to be simplified 
so that an analytical solution can be attempted. For aglsymmetr Ic , . 
isothermal. Incompressible stead, flow under weightless conditions, 
the governing equations In cylindrical coordinates can be written: 

Continuity: 


^ ' (ru) + ^ 0 

r , r •' /, 


( 1 ) 
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T 


Momentum: 


r Component : 


/ . r>ii\ ;jp J '> r I j. ^ 

»'r '' SV ■ • Vr * "te I'r Tr ^ , 2 


2 Component: 


/ . 3v\ 

'■''I?] 


5P ^ 

- 3l ^ ^ 


19 / Dv\ ^ o' 

r ar ^ 9r) 


Boundary conditions are required on the flat plate, along the 
axis of symmetry, at the nozale exit, on the free surface, and at 

r “ L. 


On plate 


u = 0 


V = 0 


on z ■ 0, all r 


Along axis of symmetry 


u = 0 


At the nozzle 


on r = 0, all z 


V » -V 


u = 0 


on 0 £ r < Rq, z = H 


At r = L 


u = u(z) 
V = v(z) 


on r = I,, 0 < z < f(L) 


On the free surface, denoted by - f(r), two boundary condi- 
tions are required since the free surface position Is an unknown to hi- 
determined as part of the solution. The details of the calculation 
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for tlK' boundary condltlonn aloiiK the froo nurface ran ho found Jii Ap- 
iH-ndlx A. (Hoc eq«. (A.IO) nd (A. 14).) 

On tlie free aurface 


2 (u^ + v^) - ~ ^ 

^ pr dr P 


V ^ (f y 


1 y2 _q 

2 ^ ~ OR 


an - f(r) 


-u - + „ . 0 


on z = f(r) 
s 


In equation f4) , the no-flow and no-slip boundary conditions are 
applied at the wall. Equation (5) is a statement involving the known 
geometrical symmetry of the problem, while equation (6) Imposes an 
initially uniform velocity profile on the incoming jet. Equation (7) 
simply states the velocity distribution as the liquid leaves the con- 
trol volume. Equation (8) is a statement of conservation of mechanical 
energy along a streamline, while equation (9) states that the normal 
velocity component on a streamline is zero. The second terms on the 
left and right sides of equation (8) ai*e the contribution of surface 
tension to the mechanical energy balance. 

The solution of the problem can be greatly facilitated by simpli- 
fying equations (1) to (8). Specifically, the method of obtaining the 
minimum parametric representation of a problem will be employed in 
order to simplify the governing equations. This method is described 
in detail by Krantz (23) and is the most systematic approach for .scal- 
ing the governing equations. The Initial stop in the minimum paramet- 
ric representation method is to form dimensionless variables hy Intro- 
ducing characteristic scale factors for .all dependent and Indopondent 
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variables. The unknown scale factors are defined as Uq, V„, r^, 
Po* ^o* 

Dlmenslanlcss variables are now defined as: 


, V* “ r* ■ a* 

U * V r 

''o '■o 


. P* - f* - 7- 
ro 


( 10 ) 


Introducing these dimensionless variables into the differential 
equations and boundary conditions, and arbitrarily making the coeffi- 
cient of one term In each differential equation and boundary condition 
equal to unity, results in: 

Continuity: 


P ^ (r*u*) = 0 
az* V r r* 3r* ^ ^ 

o o 


( 11 ) 


Momentum: 

r Component: 

^^0^0 ^ 3u* ^^0^0 ^ 3u* _ ^ ^0^0 

r p 3r* u 3z* ~ r pU 3r* 

o 00 


1° _i_ 

2 3r* 




9^v * 

3z*^ 
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z Component: 


^ u* i2il^ + V* iv* 

r V “ 3r* 3z* 

o o 
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o ?p* 
2 3z* 
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o o 
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r* 3r* 
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o' O dZ* 
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Boundary conditlonn 
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At r = L/R, 


u* = u*(z*) 


on r* = 
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V* = v*(z) 

On the free surface 
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The scale factors must now be determined. This is done by sottiiiR 
some of the resulting dimensionless groups in the ec|uatlons and bound- 
ary conditions equal to zero or unity, the groups chosen depending upon 
the physical conditions for which the equations are being scaled (23) . 
Characteristic lengths are usually determined from the dimensionless 
groups generated by the boundary conditions, while characteristic 
times, velocities, etc., are determined from dimensionless groups 
generated by the differential equations. The guidelines in deter- 
mining the unknown scale factors are: 

Do not Introduce any mathematical contradictions 
(2) Do not violate physical intuition. 

Boundary conditions. - Examining the boundary conditions, it is 
apparent that the following dimensionless groups are introduced; 


V_ 

/ 

o 


r ’ 

o 


V » r * 2 ’ r.’ 


IL 

t 

o 


r 

o 


U- 


af 


f f U f 

_o a o o c 

„2* „2 2’ 2’ 2’ z ’ V r 

V pVr r pRV o oo 

o o o o o o 


It is known that v has the range 0 to -V, r has the range 
0 to L, 2 has the range 0 to H and f(r) also has the range 
0 to H. Therefore, setting. 


* 1 


and 


( 20 ) 


Implies that 


V 


and 


o 0 

and yields two of the six unknown scale factors. 


( 21 ) 


Some of the above remaining dimensionless groups cannot bo set 
equal to one or zero without introducing contradictions. Setting 


V = V and r = f, into the above ratios, and since It would be ex- 
o o 


pected that 
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the following meaningful ratios remain: 


2 ’ 2 2 * 


Setting the second or fourth ratio equal to one or zero wo >ld violate 
physical intuition. Therefore, the ratios to be considered are 


"o’ PV^L^’ 

At this point an attempt was made to set H/z^ » 1 such that 

would equal H. This seemed logical because 0 to H was the range 

of z. However, this leads to some confusing results in terms of the 

physics. For a given flow condition, it Is argued that for a certain 

(minimum) value of H up to H *= », the flow pattern in the vicinity 

of the plate Is not expected to change. This Is shown schematically 

In Figure 2. This argument has been experimentally verified and will 

be discussed at some length In Section IV, Experimentation . The major 

point Is that H cannot be a characteristic length In the problem 

either with reference to z or f . This leaves two remaining ra- 

0 0 

tios from consideration of the boundary conditions. 


pv2l2 ^ 

Accordingly, all possible information from the boundary condi- 
tions has been obtained. Two of the six scale factors and a relation- 
ship between two others has been determined. The governing equations 
must now be examined. 
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B. Continuity Equation 

From the phyalcs of the problem, It la known that mass must be 
conserved. Hence, the continuity equation nust be valid in Its dimen- 
sionless form (eq. (11)), If the dimensionless derivatives 3v*/3z* 

and (1/r*) (3/3r*y (r*u*) are to be of the same order of magnitude. It 
Is required that 


V^r^ (23 

o o 


With V 


V and r. 


o o 
VL 


= 1 


L, it Is found that 


(24 


Solving for U^, 

o z (25 

o 

This, of course, is an equation relating two unknowns, U and z I 

o o 

is noted that the same information could have been obtained by setting 
the second of the two ratios remaining from the consideration of the 
boundary conditions equal to one. 

C. Momentum Equations 

It is the objective of this analysis to define that portion of 
the flow for which an invlscid solution is valid. For this case, the 
pressure forces are balanced by the. inertia forces. This f.ict allows 
u, to determine the scele f.ctor for the pressure, If the dleen- 

slonless pressure gradient in equation (13) is to be the same order 
of magnitude as the dimensionless inertia term. 
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^ Theae scale factors can „o» he suhstltct.d Into the governing 

Obtain the minimum parametric representation of the 
Ptoblcn. The results are aa follows, 
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( 30 ) 


(31) 


where Re is the ReyiiOlds number defined as 



The equations are now in the form in which it can be determined what 
the conditions must be such that viscous forces are not significant. 
The major parameters in this probJem are (R^/L) and Re. Since all 
the starred or dimensionless terms In equations (30) and (31) are of 
unit order, only the coefficients of the Individual terms need be 
considered in order to make statements regarding simplifications. 
Considering equation (30), it can be seen that since (R^/L) <* 1, 

both the inertia and pressure forces will be an order of magnitude 
greater than the viscous forces provided that, 

Re » 1 



From equation (31), since (R^/L)^- « 1, no new information is ob- 
tained. The governing equations will be reduced to a slmpllflod form 
of Euler’s equations of motion. 
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( 34 ) 


u* 


3v* 

3r* 
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V* 


3v* 

3z* 


3p* 

3z* 


(35) 


E. Results 

Restricting the viscous forces to be at least >.v,o orders of mag' 
nitude smaller than the inertial or pressure forces, there results 
Re > 100 


>100 (i6> 

Since (R^/L) is less than one by definition, the coefficient to con- 
sider is the second one in equation (36), since this will be the limit- 
ing one. Under the following restrictions, 

(1) Re > 100 



(37) 


the Euler's equation of motion are obtained. The equation 
Re(R^/L)^ = 100 is shown graphically in Figure 4. The line shown in 
Figure 4 separates the Inviscld from the viscous region. A physical 
understanding of the problem is made clearer by reference to this 
plot, The higher the Incoming Jet Reynolds number becomes, the 
thinner will be the boundary layer at some fixed radial position from 
the stagnation point. As seen from Figure 4, at lower values of the 
ratio R^/L (perhaps obtained by increasing the disk radius L) , 
higher Reynolds numbers are required in order to avoid viscous influ- 
ence. Finally, within the viscous region, the boundary layer will 
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grow until a radial length Is reached where It becomes equal to the 
free surface height. The effective design of an experiment Is now 
possible so that the flow can be considered essentially Inviscld, 
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IV. ZERO GRAVITY EXPERIMENTATION 
A. Apparatus and Procedure 

1. Test facility . '• The experimental investigation was conducted 
in the 2.2-second drop tower. The exact specifications of the fa- 
cility, the mode of operation, and release and recovery systems are 
described in detail in Appendix B. The drop tower provides us with 

a 2.2 second weightless environment in which to conduct the tests. 

2. Experiment . - The experiment package used to obtain the data 
for this study is shown in Figure 5. It consists of an aluminum 
frame in which were mounted the jet reservoir, disk, a 16 millimeter 
high speed motion picture camera, supply tank, backlighting scheme, 
and batteries. The major functions were controlled by on-board se- 
quence timers. 

A diagram indicating the manner in which the flow system operates 
is shown in Figure 6. This is a pressure controlled system in which 
the flow was initiated by opening the solenoid valve. Prior to the 
drop, liquid was contained in the line between the liquid supply tank 
and the jet reservoir. In addition, the jet reservoir was completely 
filled with the test liquid. 

Initially, the object was to employ two-dimensional jets during 
the experiments. Schach (38), was the only author who investigated 
two-dimensional jets, although it is impossible to determine how rec- 
tangular his jets really were. Several attempts were made to fabricate 
a two-dimensional nojizle (slot) for use in the zero gravity studies. 
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Initially, the slots built were quite crude; a rectangular hole cut In 
a section of plexiglass; a balsa wood jet sealed with epoxy cement. 
Later, slots were accurately designed in order to achieve the desired 
flow. The slot length was arbitrarily chosen as 5 centimeters and Its 
width as 0.25 centimeters. Approximately 45* tapers were made to thi> 
opening along both the . Ide and narrow sides. The results of the zero 
gravity testing were as follows: The jet contracted In the long direc- 

tion '5 cm) and expanded In the narrow direction (0.25 cm). A tendency 
to become cylindrical probably due to the effects of surface tension, 
was observed. A redesigned version of the slot was tested, which em- 
ployed absolutely no taper in the narrow direction. This also failed 
to yield a rectangular jet. When these approximate slot jets Impinged 
on various flat plates, the jet diverged radially. In ether words, it 
was Impossible to prevent spreading in the lateral direction for an un- 
constrained surface. Various methods were tried to eliminate this 
lateral spreading and to force the impingement flow to be completely 
two-dimensional. The flow was Impinged on rectangular plates having 
the same width as that of the slot, 5 centimeters. The result of these 
tests was that the liquid simply fell off or went around the plate. 

A large rectangular plate was used and the region greater in width than 
5 centimeters was sprayed with fluorocarbon since distilled water does 
not wet a fluorocarbon surface. This also failed. All attempts at 
using a two-dimensional slot jet were subsequently abandoned and axl- 
symmetrlc jets were pursued. 

A schematic drawing of the Jet reservoir, which was fabricated 
out of acrylic plastic, is shown in Figures 7(a) and (b) . The critical 
feature of the jet reservoir is the 45“ taper to a circular hole of 
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diameter D. The range of diameters studied was from 0.5 to 1,5 centi- 
meters. In addition, the transition between the circular hole and con- 
ical taper was rounded smooth. The taper prohibits boundary layer 
buildup and allows the liquid jet to exit from the reservoir with a 
nearly uniform velocity profile. Based on the experimental results of 
reference 24, a total angle of 90" is more than sufficient to insure a 
uniform exiting velocity profile over the range of Reynolds numbers 
studied. 

Sharp-edge disks, also fabricated from acrylic plastic, were 
mounted above the jet reservoir by means of a threaded rod such that 
the flat surface of the disk was at right angles to the impinging li- 
quid Jet. Sharp-edged disks were employed in the study rather than 
finite thickness disks in order to develop an analytical model for the 
flow. As was learned later, the thickness of the disk edge is only of 
importance for those flows dominated by the effects of surface ten- 
sion. The diameters, 2L, of the disks were 2.0 and 3.0 centimeters. 

A schematic drawing of the disks is shown in Figure 8. 

3. Test liquids . - Two test liquids were employed, anhydrous 
ethanol and trlchlorotrifluoroethane. Their properties at 20® C are 
listed in Table 1. No attempt was made to correct the fluid properties 
for temperature changes. It is noted that both of these test fluids 
possess a nearly 0“ contact angle on an acrylic plastic surface. How- 
ever, this was not the reason why they were chosen as teot fluids. 

They were chosen because of their relatively low viscosity and avail- 
ability. 

4. Test pr oc edure . - Prior to a test run, the jet reservoir, disk 
and supply tanks were cleaned ultrasonlcally with a mfld detergent. 
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After these parts were rinsed with methanol, they were dried in a warm 
air dryer. The supply tank was subsequently filled with the test liq- 
uid and the Jet reservoir was filled by pressurl^.lnB the supply tank. 

This procedure eliminated air bubbles from the lines ensurinR accurate 
flow rates. After the jet reservoir was completely full, the supply 
tank was sealed and two accumulator bottles (not shown in Fig. 5) were 
pressurized with gaseous nitrogen to a predetermined value. The ac- 
cumulator bottles were designed to be of such a volume that no appre- 
ciable pressure drop occurred during the drop. 

Electrical timers on the experiment package were set to control 
the initiation and duration of all functions programmed during the 
drop. The experiment package was then balanced and positioned within 
the prebalanced drag shield. The wire support was attached to the ex- 
periment package through an access hole in the shield (see Fig. B3(a) 
in Appendix B) . Properly sized spikes tips were Installed on the drag 
shield. Then the drag shield, with the experiment package inside, was 
hoisted to the predrop position at the top of the facility (Fig. Bl) 
and connected to an external electrical power source. The wire support 
was attached to the release system, and the entire assembly was sus- 
pended from the wire. After final electrical checks were made and the 
experiment package was switched to internal power, the system was re- 
leased. After completion of the test, the experiment package and drag 
shield were returned to the preparation area. 

B. Experimental Results 

X. G eneral considerations . - In addition to the measurement and 
observation of steady state liquid flow patterns, two separate phenomena 
were observed. First, no circular hydraulic jump occurred during any 
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of the teats even though they are a common occurrence under normal 
gravity conditions. Secondly, the initial Impact of the jet upon the 
solid surface provided another unusual phenomena, that of tli^ rebound- 
Ing liquid droplet. At high flow rates, the jet broke up prioi to im- 
pinging upon Che disk. The fiist droplet tended to impinge upon and 
stick to the disk, spreading as it did. However, the second droplet 
impinged and rebounded off this wetted surface sometimes into the in- 
coming liquid jet. This provided no serious problem with the attain- 
ment of a steady state flow pattern since this all occurred during 

the transient phase. 

The jet generally appeared to go through three phases during the 
Impingement process. The Initial phase. Including the droplet pinch- 
off and subsequent impingement, was termed the transient phase. After 
a certain period of time, a steady state flow pattern was achieved 
from which the free surface shapes were measured and observations 
were made. A third phase was reached shortly after the jet had 
reached its equilibrium configuration. The flow pattern developed an 
instability. Initially, the instability started from the jet which 
began to oscillate, and then spread to the plume. Since the time over 
which the jet and flow pattern appeared stable and smooth was finite, 
steady state data could be cbtained. It was observed that the time 
before breakdown was inversely proportional to the back pressure dur- 
ing the flow. At 1 psia, for example, the jet remained stable for 
1.5 seconds while at 10 psi, It was stable for approxim.ntcly 0.4 sec- 
onds. An attempt was made to lengthen the time over which stability 
occurred by packing the nozzle chamber with steel wool. This appeared 
to be effective in improving the overall stability but had the negative 
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effect of IntroduclnR a low lovol |H>rtiu hal Ion thrmij'lidui I he teal aiitl, 
thuH, wa« not employed dnrlnp, any f)f the f'xpnrimenitfi , 

2* of. ln*ti»jht. ~ Si'ver/il I <•»! s werr* coiKhici e<l Ini- 

tially to (leteniilne the effect of il, the distance hetween I lie iin/.zle 
and the disk, on the experimental flow pattern. i| was delennliied tliat 
there ia no effect on the liquid flow provided ihal the ratio of the 
nozzle height H to the jot diameter is great er than 1. As a ri'sult, 
all zero gravity experlment.s were conducted in order to eliminate this 
effect. This fact, which was determined experimentally, supports the 
argument made In Section III, Order of Magnitude Analysis , concerning 
the assertion that H could not be a scale factor for the axial co- 
ordinate. 

3. Steady state flow patterns . - The approximate steady state 
flows for three different Jet velocities are shown photographically In 
Figure 9. The direction of flow of the liquid jet Is vertically up- 
wards. The threaded rod and bolt observed in the film clips Is the 
disk holder which connects the sharp-edged disk assembly to the rig 
frame. The jet velocity increases from left to right in the figure. 
Three distinct classifications of flow patterns were observed to occur 
and are shown labeled in Figure 9 as surface tension flow, transition 
flow, and inertia flow. Surface tension flow (Fig. 9(a)) is defined 
as that flow in which the liquid flows completely around the disk with 
no separation occurring from the disk edge. In transition flow 
(Fig. 9(h)), surface tenslua and inertia forces are botli important. 
Transition flow is defined as flow in which separntion occurs from tlie 
disk and the resulting liquid sheet either eollecl;. upon itself form- 
ing an envelope or has the tendency to do so. Inertia flow (I-’ig. 9(c)) 
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is doilnofl ns tlint flow In whlcli the llcpild ,'ie|»ai /jtos frfini I lie dfr.lt 
with no liquid turniii}^ towm dti the jid; eonterUiu* ntf oinpt 1u(' to fori;i 
an onvi'lopc. The flow p/ictern nhown for transition flow Is noi rrally 
Iho fitn.’idy ntato flow pattern one woul<l t'xp<>ct if fiti'.idy stale eond I 
lions could have b»'en readied. T)k> reason for this is .is follows; 

Some liquid Is .always travtsUii}' toward tlic disk from llie polnl .ac 
which the envelope meets. This licjuid flow strikes tlur b.ark of the 
baffle and subsequently disrupts tlie Initially formed rmvelope. Tlie 
recirculation flow then is a strong function of the geometry of the 
disk holder, an uncontrollable parameter. This important distinction 
means that an analysis for free surface shapes would li.ave to account 
for the geometry of the disk holder in the transition region. The 
surface tension flow is generally slow and, as a result, does not 
quite reach a steady state configuration on the back side of the 
disk. The inertia flow tests (represented by Fig. 9(c)) always 
reached steady state. 

The experimental tests were conducted in the inviscld region of 
Figure 4. Depending on the particular ratio of (R^/L) , there exists 
a minimum Reynolds number below which the runs can no longer be con- 
sidered as viscous-free. The experimental results are listed In 
tabular form In Table 2 in which all the Important parameters as well 
as the flow classifications are contained. One additional parameter, 
the Weber number, Is listed in Table 1. As will be shown in the next 
■ ection, when the Reynolds euniber is no longer a parameter to con- 
sider for flow classifications, tiie Weber number and ratio /h) re- 
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main. The ' her number pV i« h.isically the ratio of inertia 

to surface tension forces. In the flow category column, S Indicates 
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surface tension flow, T Indicates transition flow, and I is inertia 
flow. Finally, it is noted that the designated flow classification 
for some cases, pattlculariy those bordering transition or inertia 
flow, could easily fit Into either category. 

fiya . vfty result s. - The data contained In TabJe 2 Is 

shown graphically in Figure 10. The lines indicated in the figure 

were faired in by hand and separate the various flow classifications. 

It is observed that at any particular value of the ratio (R^/L) , the 

flow classification is dependent only on the system Weber number. 

Two critical Weber numbers occur at a constant value of (R^/L). The 

lowest critical Weber number separates the surface tension flow from 

the transition flow while the higher critical Weber number separates 

the transition flow from the Inertia flow. In addition, the critical 

Weber number between regimes was found to decrease as (R /l) was in- 
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creased. 
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V. POTKNTIAI. F()RMtJ;,AT10N 

A. Covornlng Equations and Boundary Conditions in Physical 

Plane Including Surface Tension 

In Section HI, Order of Magnitude Analysis , it was shown that 

at any particular value of R /L if Re > Re the flow in the re- 

o cr 

glee of the disk can be considered as viscous free. It is further 
assumed that the Jet will continue to remain viscous free after leav- 

i 

Ing the disk. There will be no shear stress between the exiting radial 

! 

jet and the ambient air surrounding it. i 

i 

1- General formulation . - Consider the flow of a circular liquid 
jet impinging normally on a circular disk as shown in Figure 11. R is 

O i 

the jet radius, I. is the disk radius, and H is the distance between ' 

the jet reservoir and the disk. The incoming jet velocity Is given as 
V and the initial velocity profile is assumed to be uniform. There 
will be two free surfaces involved. The upper free surface is defined 
z = f. (r) and the lower surface as z = f^(r). In addition, a 
third surface is required for the complete formulation of the boundary 
value problem. Initially, this surface was chosen as the straight lino 
FE shown in Figure 11. However, this .-roved to be inconvenient since 
FE Is arbitrary and, thus, has no known boundary condition. it was 
found convenient to instead choose the surface z = fj(r) to he ;i sur- 
face along whicli the velocity potential is constant. Various points 
in th(’ physical plane li.ive bec-n designated with Ic-ttevs riiy.inr, from 
A to C for ctuiven i t'n<H- , In cyllndric.il coordlnal e.s, the govetniiU’- 
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equations and boundary conditions in terms of primary variables (u,v) 
are given as follows: 

Continuity: 

1 A (,„) + |2 . 0 (3 

r 3r ' 3z 

Momentum: 

r Component; 

/ 3u . 3u\ 3p (3 

z Component: 

( 3v . 3v’\ _ ^ 

The following boundary conditions are applied, 

On DC 

v " 0 On 7. = 0,0<r_<L 

On AB 

u“0 Oli“£Rg, 2 *H (' 


On 

BC 



3v 

3r 

= 0 

H 

II 

O 

w 

0 < z 

On 

AG 



(u^ + 

v^) 

<3 d ( 


- pr dr r 

^1 + f 


1 

2 pR 


on z = f , (r) (44) 

S 1 


-uf + v “ 0 
On DE 


on z “ f, (r) 
s i 


1 . 2 . 2 . 

2 (u + V ) 


_L Jl 

pr dr 


1 

2 pR, 


on X r,(r) 

s - 
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-uf2 + V “ 0 
On GE 

u + vf^Cr) * 0 


on Zg “ 


on 2 " ^^(r) 


The derivation for the boundary condition along GE is as fol- 
lows. /.Irng GE, the velocity potential ^ is constant. By defini- 
tion, the velocity vector ^ is normal to an equipotential surface. 


This means that 
^ X n ■ 0 


on z = fn(r) 


The unit normal to f^Cr) is 

. - 

VJFTT 

where i is the unit vector in the radial direction and j is the 
unit vector in the axial direction, and 
7 “ ui + vj 

Application of equation (49) results In 
u + vf '(r) 

— ^ a 0 

yfp + 1 

From which equation (48) follows. 

2. Introduction of Stokes Stream Function . - The primary variables 
contained in the governing equations are the scaler velocity components 
u and V. The fact that two functions are required to describe one 
vector field is cumbersome. As shown in the theory of hydrodynamics, 
the number of functions can be reduced for several important cases, 
one of these being axlsynmetric flow. A function ip, defined as 
Stokes stream function, can be introduced which automatically satisfies 
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the continuity equation. With the additional requirement of Irrota- 
tlonality, the governing equation in terms of Stokes stream function 
will result. According to Chan (7), Stokes stream function is a mathe- 
matical device used to describe the flow and has the following proper- 
ties. First, when the stream function is set equal to a constant, it 
results in different annular stream surfaces in axisymmetrlc flow. 
Secondly, when it is differentiated properly, it yields the velocity 
components. Thirdly, taking the difference between the values at two 
adjacent stream surfaces yields the flow rate. 

Starting with equation (38), the continuity equation for axisym- 
metric flow, a guess is made at what u and v are in order to sat- 
isfy continuity identically. Assume 

(53) 

r 3z 


And 


(54) 

r 3r 

Substitution of equations (53) and (54) into equation (38) shows that 
continuity is identically satisfied by these two guesL.es. In addition, 
the assumption is made that the flow is also irrotational . As a re- 


sult 


Curl V “ 0 

For axisymmetrlc flow, this can be written, 

. lx = 0 

3r 

Replacing u and v in the above by thei’- relat im-shiiis to 
in cVie governing equation for axisymmetrlc flow in i orms ol 
■stream t unction 


(55) 


(5ft) 


ij) results 
S l okos 
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. i a + rt . 0 


3* Deriva tion of boundary conditions In terms of it/ . - Since BC, 
CD, and DE are all a part of the same btreamline, they must all have 
the same value for Stokes stream function. Let us arbitrarily set 
that value equal to zero. The value of alon • AB and AG can be 
calculated from equation (54) . 

On AB 


-V = i M 

r 3r 

Integrating this yields 



On AG since AG Is a line of constant 



(58) 


(59) 


(60) 


In addition, equation (44) applies. Substitution of equation (53) 
and (54) into (44) yields 



(62) 


It is noted that the second boundary conditions on the free sur- 
faces, namely + v 0 and -uf^ + v = 0, are fully equivalent 

to the specification of the value of Stokes stream function along that 
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surface. Since ^ » i|;(r, 2 ), It can be expanded and dip sot equal to 
zero along AG and ED. 

On GE, equation (48) becomes 

“ **■ If "0 on z - f 3 (r) (63) 

Npndlmenslonallzatlon of governing equations and boundary con- 
ditions. In physical plane. - The governing equation In terms of Stokes 
stream function (eq. (57)) and the boundary conditions (eqa. (58) 
to (63)) are now put into dimensionless form by introducing arbitrary 
scale factors. Let the scale factor for the stream function be 

o 

Let dimensionless quantities be represented by stars, l.e., ip* Is 
dimensionless. The results of this manipulation are shown In Fig- 
ure 12. Three parameters appear in the specification of boundary 
conditions. They include the Weber number, We, and the dimensionless 
length ratios L/R^ and H/R^. Recalling the arguments in the Order 
Magnitude Analysis section, K/R^ is not really a parameter pro- 
vided it is larger than some minimum value. The dimensionless velocity 
components can be calculated from the following expressions 


_L 

r* 3z* 

(64) 

J;. lit!* 


r* 9r* 

(65) 


The procedure for solution of the boundary value problem as set 
up in dimensionless form in the physical plane (Fig. 12) would be as 
follows: Initially, realistic variations for f#(r), f*(r), and f*(r) 

are assumed. Using only one of the two given boundary conditions on 
AG and GE (i.e,, ip* • 1/2 and ip* » 0) solve for i/»* using finite 
difference methods. With the initial solution for i/* check the va- 
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lidlty of tho flecond of the two boundary conditions on AfJ and OH. 

If the boundary condltlonn are not satisfied, new variations In f*, 
f*, and must be assumed, A second iteration to ip* must be ob- 

tained and so on. One of the serious drawbacks of this outlined it- 
eration scheme is the lack of knowledge concerning how to update as- 
sumed values of fj, fj. and based on previous solutions. In 

other words, there is no logical way in which to make changes to the 
shape of the initially assumed control volume. 

5, Surface tension dominated model . - For Weber numbers between 
5 and 30 (depending on the ratio Rq/I.) , experimental data shows that 
the resulting steady-state flow pattern is surface tension dominated 
(see Fig. 9(a)). By previous definition, surface tension flow is de- 
fined as that flow in which the liquid flows completely around the 
disk with no separation from the edge. It is the intent to model this 
flow in order to solve for the theoretical free surface shapes and ve- 
locity profiles. Assuming axlsymmetry, the physical plane model is 
shown in Figure 13. In the model, at some cross-section far down- 
stream, the flow is assumed to approach the Initially uniform flow it 
possessed at AB. The exiting plane is denoted by GE in the model. 

C and C are both located at r=0.z*0. C is located on top 
of the plate while C is on the bottom. The free surface fj^(r), 
is not assumed to possess mirror symmetry about the z = 0 position. 
In cylindrical coordinates, the governing equations and boundary con- 
ditions in terms of the primary variables (u,v) are given as follows: 
Continuity: 

(ru) + 1^ » 0 


1 J_ 
r 9r 


( 66 ) 
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Momentum; 

r Component 
/ 3u . 3u\ ^ 

“ ^ ij) ' • at 

z Component 
/ 3v , 3v\ ^ 

» ^ ' al) ' ■ at 

The following boundary conditions are applied 
On DC 
V - 0 
On DC' 

v”0 on z = 0,0£r£L 

On AB 
u “ 0 
On BC 


3r 
On GE 

u •» 0 

On AG 


on z“0, 0<r<L 


0^r<^R^, z = H 


1^*0 on r»0, 0<z<H 
3r - 

On C'E 


-^»0 on r-0, -H<z<0 
3r ~ 


0 < r < R^, z *= -H 


1.2.2. ad 

2 



2 pR„ 


(67) 

( 68 ) 

(69) 

(70) 

(71) 

(72) 

(73) 

(74) 

on z = f (r) (75) 

S X 


and 

-uf^ + V = 0 on Zg = ^1^^^ 

Equations (75) and (76) represent two distinct pieces of information 
concerning z^ = fj^(r). It is noted that fj^ can be eliminated from 
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the set of equations by means of simple substitution. Tills fact will 
have Implications later in the development leading to a major simpli- 
fication. The following expression would result, 



Direct substitution of equations (53) and (54), the relations between 

the velocity components and Stokes stream function, Into the governing 

equation and boundary conditions for the surface tension model results 

in the formulation shown in Figure 14 after nondlmenslonallzatlon. 

Similar to the general formulation in the last section, the various 

lengths in the problem are scaled with respect to R and the scale 

o 

2 

factor for Stokes stream function is -VR . 

o 

B. Inverse Plane Formulations 

The procedure for solving for the free surface shapes in the phys- 
ical plane has been outlined previously, Tne difficulties encountered 
when making adjustments to the free surfaces between iterations and the 
lack of a logical manner in which to make the adjustments have been 
cited. This would be a time consuming task even in the absence of sur- 
face tension. A computerized scheme is sought which offers the possi- 
hllity of achieving the free-surface results with a minimum of computer 
iteration time and user interaction. 

1. Transformation formulae . - An alternate approach to the phys- 
ical plane solution is discussed In detail by Jeppson (19), Tn his 
article, Jeppson discusses a transformation technique into what Is de- 
fined as the inverse plane. The coordinate axes in the Inverse plane 
are the velocity potential and Stokes stream Tunction J'. The .'id- 



39 


vantage to using the (Juji apace, or inverse plane, is that the free 
surface lies along a line of constant and Is, therefore, at a 
known position. Of course, one must pay the price for knowing the 
position of the free surface. As will be shown shortly, the governing 
equation Is no longer linear as It was In the physical plane. 

Intxoductlon of ve locity potential . - The use of a scaler func- 
tion defined as the velocity potential has been used previously In the 
>tpeciflcation of the boundary GE (the exiting plane). The continuity 
equation for steady incompressible axisymmetric flow Is given by equa- 
tion (38). The flow is also assumed to be Irrotatlonal such that the 
condition given by equation (56) Is also valid. Equation (56) implies 
that there exists a scaler potential function cj such that 

V = grad * ( 73 ) 

front which it follows 

u = 

9r (79) 

and 


94 

V St - I . 

9z 

Substitution of equations (79) and (80) 
for steady axisymmetric flow In terms of 


(80) 


yields the governing equation 
the velocity potential 



1 91 9 ^ 

9r2 


0 


(81) 


Rs^Atlonship betwee n 1 and 4 . - It can now bo stated that the 
relationships between the velocity potential and Stokes's stream func- 
tion are given as; 


u 


= _ A 11 

9r r 9z 


and 


( 82 ) 
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V = H . i .Silf 

r :)r (H'l) 

- In order to formulate the prohJem 
In the physical plane, r and z miiat be known a« functions of ,/. 

■uul 4 >. In other worda, we are attempting to reverse the roles of the 
dependent and independent variables in the Inverse plane. The required 

relationships are obtained noting that If ^ = Kr,z) and 

then there exist inverne functions r » r(^^) and z = as shown 

in reference 13 such that 


Or _ 

= _ i 



3(fi 

J 

Oz 

(84) 

_ 1 . 



OiP 

J 

Or 

(85) 

Oz ^ 1 3ip 



9^ 

J Or 


(86) 

Or ^ 

. i ii 



Oip 

J Oz 


(87) 

where 

the Jacobian J Is defined as 


3 ji 




Or 

Oz 

■ ) - (HXIf) 

J « 

3^ 

Oip 


Or 

Oz 



Fro» the above set of e,ueclo„e follow two Itiportont relatlooe. Sub- 
atltutlng into equation (82) the value, uf H/)r and 3*/3a obtained 
from equations (84) and (85) yields 
Bz 


dip 


i lE 

r 3(|> 


(89) 


sinllerly, eubatltutlon of the values of s,/ 3 t „„d 

tlons (86) and (87) Into equation (83) yields 

Or ^ 1 

Oip * r O ir 


( 90 ) 


I 


For Bubsequont rolntloriH, It ia necoHsary to PxprcBH tho .lar.obJnn 
In torma of the inverae function r(ij»,(|-) mul Suhatlt»ttinH 

equatlona (84) to (87) into equation (88) reanlta in 


2 17 \ . pi Y-^' 


.1 = J 


Now, using equations (89) and (90) in the above yields. 


5. Governing equations for inverse functions . - Differentiating 
equation (89) with respect to \J/ yields 

liz = i ^ J: M (93) 

“ r 9i|; 94. J.2 9i|) 9* 

and combining this with the derivatives of equation (90) with respect 
to 4* which is 

1 + i A (94) 

it is found that 

A _Lifz _l,9r^ 

3*2 ,3 9* H - ,2 3^2 ^2 3* 3* 

By using equations (89) and (90), the terms involving derivatives 
in equation (95) can be expressed entirely in terms of z, giving tl)e 
important equation, 

r2 . .2 

9(ti^ 9(j)^ 

On the other hand, differentiating equation (89) with respect to 
i and equation (90) with respect to ij) and «;ombinlng the result leads 
to the eqtiation 


( 


^ ^ J . AJ. 1 / ^ _ .1 Ar ■')?; 


Hy subHtltuttnn f..i frnm o<,,i.u I,,,, ( 90 ), (Ik- I ol] owi o)’ cim, I Jon 

for r ( <fi I ) 1 /) oj) I: n :l f md J 

+ ± A^A . I f^r y . J ^ 

For further dl.sniHnJon, cquiUJon (9f,) wijj |„, rM errt d to nr t he 
a equation while equation (98) wJ,U he railed the r equation. m,th 
of these equations are nonlinear with the nonlinearity in the z equa- 
tion also involving r. It is further noted tl.at both of these equa- 
tions are of the elliptic type, as shown in reference 19. This means 
that boundary conditions are required for all boundaries In the flow 
region. Since the r equation (eq. (98)) only Involves that variable, 
a solution to the problem begins with its solution. 

One final point Involves the fact that the variables discussed in 
the above Inverse transformations are all dimensional. It is noted 
that nondimensionallzation of the above equations similar to the nen- 

dlmenslonallzation done in the physical plane allows us to recover the 
same equation. 


r* = 


. 1 
o 


where z 


where r 


wliero i/) 


“ ■■ t 


where <i = 
o 


i 


I 
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The reaulCH of this substitution yields exactly the same governlnK 
equations and transformation formulas in starred notation. In other 
words. 


r* 


0 


3^z* . 3^z* 
2 ^ 2 


^ 3z* 3z* 
^ 3i|i* 3(j)* 


and 


(100a) 


3^r« _1_ a^r* 1 /8r*\^ X /3r*\^ = 0 (100b) 

3(,*2 r*2 3^.2 ■ r.3 

6. Infinite flat plate excluding surface tension . - The first 
problem to be examined in which we employ the inverse transformation is 
the flow of a circular liquid jet normal to a flat plate. At this 
point, we drop the starred notation to designate dimensionless quanti- 
ties. It is assumed that all quantities appearing henceforth are di- 
mensionless unless otherwise stated (l.e., r, z, (|), are now dimen- 
sionless). For the case in which surface tension is excluded, the 
physical plane model Jeplcted in Figure 12 must necessarily change. 

The model in Figure 12 is for the general case of a finite disk in 
which surface tension effects are Important. For the case of an in- 
finite flat plate, since the free surface ED no longer exists, the 
designation D does not appear in the model. Furthermore, for the 
case in which surface tension forces are unimportant, the boundary con- 
dition along the upper surface (referring to Fig. 12) can be written 
as follows: 


On AG 



( 101 ) 


This expression is obtained by allowing the Weber number to approacti 


T 


T 


T 
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infinity. The reasoning for letting the Weber number become large la 
because it is the ratio of inertial to surface tension forces. The 
<'ase we are examining is one in which the surface tension forces become 
vanishingly small or the Inertial forces becoming large. Finally, the 
boundary GE will be nearly vertical implying that the velocity there 
is purely radial. 

7. Derivation of boundary conditions in inverse plane for r 
formulation . - 
On AB 


1 2 

In the physical plane, tp = -^ r 

In the inverse plane, r = -/lip 
On BC 


( 102 ) 
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In the Inverse plane, using equations (82) and (83) , this becomes 



and upon using relations (85) and (87) 

[(It/ ^ ml = ^ 

Substituting In for from equation (92) and then using (89) 

and (90) yields the. final form 



The results of the above calculations are shown in Figure 15(a) in the 
Inverse plane. For the case of flow normal to an infinite flat plate, 
the r values can now be completely found once we solve the nonlinear 
governing equation. The variable z does not appear anywhere In the 
formiilation. In addition, the boundary conditions for the z formula- 
tion were also derived and are indicated in Figure 15(b). It follows 
that the z solution canv ot be obtained until after r values are 
known since the r values are required on the free surface boundary 
and within the Interior flow. The boundary conditions for the z 
formulation were derived as follows: 

8. Derivation of boundary conditions in Inverse plane for z 
formulation. - 


On AB 

In the physical plane, u <= 0 which implies ° ^ 

From equation (85), " 0 

This Implies that z » z(0) alone and since 

AB Is a line of constant 

in the inverse plane z = constant 


( 107 ) 


r 


T 
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On BC 

In physical plane, u « 0 


32 


From which it follows, in the inverse plane “ 0 


(108) 


On EC 

In physical plane i{i 
In inverse plane z 
On GE 

In physical plane, v *» 0 


= 0 ] 
= 0 J 


(109) 




From equation (83) , it follows that “ 0 


} ( 110 ) 


3z 


and from equation (86) , we have In the inverse plane ® 0 


On AG 


In the physical plane 


i M • («)■] 


In the Inverse plane we can take over the equivalent ex- 
pression given by equation (106) ^ 

Using equations (89) and (90) to eliminate the derivative 

2 v2 


> ( 111 ) 


involving r results in r‘ 


i/ii) + 

\ 34 '/ / 


= 1 


9. Finite plate excluding surface tension . - A second problem to be 
formulated is concerned with the impingement of a circular liquid jet 
normal to a disk of finite width. For the case in which surface tension 
effects are neglected, the model depicted in Figure 12 will vary only 
slightly. The only difference being that the boundary conditions along 
the top and bottom free surfaces are now expressed as 


I 


I 
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? [(It/ ^(11)1 


The r iormulatlon In the inverse plane Is shown In Figure 16(a). The 
only boundary condition that bears some explanation is that along GE, 
the exiting plane in the physical plane, from Figure 12, 


3r ^3 ° 


Using equations (84) and (86) this can be written 
0 


3<j> 3(|> ^3 


Finally, using equation (90) to eliminate 3z/3^, 


3r 3r ^ 

9<|> 3i/> ^3 “ ® 


( 112 ) 


(113) 


(114) 


or 


3r 

3 ^ 


The z formulation in the inverse plane is shown in Figure 16(b). 
Here again, only the boundary condition along GE bears some explana- 
tion. Continuing with equation (115) and using equations (89) and (90) 
we obtain 

<n) - (ifXf ) ■ » (U6, 

The formulation of the problem of flow of a circular liquid jet 
normal to a finite disk In the inverse plane is now complete. The r 
solution (Fig, 16(a)) can no longer be obtained independent of the z 
solution due to the exiting jet boundary condition. Thus, the problem 
will necessitate a simultaneous solution of two partial differential 
equations of the nonlinear elliptic type. 

PTate including surface tension . - A third problem that 
is formulated but not solved due to the complexity of applying It to n 
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real system involves the normal Impingement to a finite disk in which 
both surface tension and inertial forces are important i These would 
Include all the flows experimentally labeled as transition flows . This 
represents a complete inverse formulation for the exact physical model 
represented in Figure 12 • The results of the inverse formulation for 
the r solution are shown in Figure 17. The way in which the inverse 
boundary conditions along the free surface were derived will be made 
more clear In the next section. 

11. Surface tension dominated model - finite plate . - The final 
problem to be examined involves the surface tension dominated flow, 
described in Figure 14. The Inverse plane formulation for the r so- 
lution is indicated in Figure 18(a). The boundary conditions are the 
same as derived previously with the exception of the free surface 
boundary. Along AG in the physical plane, 

(117) 


(118) 


(119) 


JL f/'af + /Ilf 1 . A i -1 

2^2 [1^32/ [dr/ J We r dr 




+ 1 


1 _1_ 
2 “ We 


In the derivation of equation (106), it was shown that 

(s?) ' Isf j 

In addition, using equations (84) and (86) , 



3z \ ^ X l-M 



2 Dr 


Therefore, the inverse boundary condition along AG becomes 


I 


I 


1 




J 1 « 
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1 

2 


m ^ 'HU) 

which can be written 


lid 

r 

We r dr 

J 1 + 1 


V 2 \3r/aW 


Lf r J 

r 

1 


1 . 

We‘ 


1 

2 


-L i -i- 

We r dr 


(H) ^ 'Hit) 

Multiplying through by 2, and rearranging yields, 


J/armf + 

\V3r/3'l'/ 


1 - -L 

2 Vie 


( 120 ) 


( 121 ) 



1 


f 

2 2 1 d 

2 

r 

•N 

y " We ^ We r dr 
< 


> 


( 122 ) 


There has been no sytnmetry assumed for the r formulation, only 
axisymmetry. It is noted that the r formulation does not Involve the 
variable t. This means that the r and z solutions can be obtained 

independently. The z formulation is Indicated in Figure 18(b) . 
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VI. CENTRAL FINITE DIFFERENCE REPRESENTATION 
The finite difference operators for the nonlinear r and z 
elliptic partial differential equations resulting from the inversion 
were put into difference form. It was found from experience that 
considerable flexibility resulted if the difference equations were 
derived employing rectangular mesh. The reason for this becomes 
clearer as we progress into the numerical solution. Suffice to say 
that thii allowed us to control the size of the flow region. Square 
meshes were originally attempted, but only led to solutions in the 
cases where the mesh sizes became vanishingly small. 

A. Formulation 

1. npneral considerations . - The partial derivatives appearing 
in the r and z governing equations are replaced by algebraic cen 
tral finite difference operators. A complete derivation of the var- 
ious operators is contained in reference 35. The notation used in 
this report is that shown in Figure 19. If the finite difference 
analogy of 3/9r is 6/6r, the r derivatives can be written at 

point 0 as follows, 



ri - r^ 

(123) 

64 * 

2 "A(^ 


6r 

^2 " 

(124) 


2 



+ r^ - 2r^ 

025) 

64>^ 

A(Ji^ 
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Similarly, the central finite difference ropreaentatlona for the 7. 


derivatives are, 

fiz « ^1 ~ ^3 

6<{. ” 2 Ai|> 


( 127 : 


iz ^ ^2 ~ ^4 

6il» 2 At|/ 

^,2 



z- + z, - 2z 
2 4 o 


(128; 

(129: 

(i3o: 


2. Interior nodal points . - If equations (123) to (126) are sub- 
stituted into the r governing equation (eq. (98)) with a - 
the finite difference expression for all interior nodal points is ob- 
tained, 


< - f <-^2 * + 'o 


— - — (r - r )^ 
2 8 ^ 2 
a 


~ , 2 ^*^1 "** ^^3^ „ 2 ^’^l ~ ^3^ ” ® 

Za ooi 

In addition, if equations (127) to (130) are substituted into the z 
governing equation (eq. (96)), the finite difference expression for all 
interior nodal points is obtained 


J?7V\ * 2/rV Vi) 


+ z^) 


I 
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where, es In the derivation of equation (131), 


a 


M 

Ai{( 


(133) 


B. Excluding Surface Tension 

1. Infinite flat plate . - The finite difference representation 
^or the infinite flat plate Is shown in Figures 20(a) and (b) . The 
algebraic expressions for the boundaries la derived by simultaneous 
application of the governing equation and boundary conditions at a 
fixed point. This application Involves a fictitious point, f, outside 
the boundary which is subsequently eliminated. Point G represents a 
special point In the finite difference representation for the z solu** 
tlon since It Is a part of two separate boundaries. The governing z 
equation (eq. (132)) was applied at point G which resulted in two fic- 
titious points. Then the boundary conditions along both AG and GE 
were applied at point G. This allowed the elimination of the two fic- 
titious points from the resulting finite difference expression. De- 
tailed calculations for the boundaries are contained In Appendix C. 

2. Finite plate . - The Inverse formulation for the finite plate 
problem (Fig. 16), Is shown In difference form In Figure 21. For this 
r.^se, the difference operator along GE Is more complicated than in 
the Infinite plate case. In fact, both G and E represent special 
points in the formulation. However, one of these, point E, Is speci- 
fy as a known position (r = constant). At point G the equation to 
be satisfied in the r formulation is shown at the top of Figure 21(a). 
In the formulation (Fig. 21(b), both points G and E arc special 
points. The following equations hold there. 
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At point G 



Detailed calculations for all the additional boundaries encountered for 
the finite plate case .';re contained in Appendix C. 

C. Surface Tension 


The difference formulation corresponding to the Surface Tension 
model (shown in Fig. 18) is indicated in Figure 22(a) and (b) . The 
only boundary condition that must be explained is the one along the 
free surface AG. Details of this calculation are contained in Ap~ 
pendix C. It is noted that (see Fig. 22(a)) the r solution can be 
obtained independent of the z solution. 




VII. DISCUSSION OF NUMERICAL TECHNIQUES 

Since the govenilng equations for the r and z formulations 
are nonlinear In the Inverse plane. In general, the finite difference 
operations at the Interior and boundary points will be nonlinear (sec 
eqa. (131) and (132)). For the Infinite flat plate case, as described 
In Figure 20, the solution begins with r = r(ij),4)) from Figure 20(a). 
Secondly, with a knowledge of the r solution, the z formulation, 
shown in Figure 20(b), Is solved for z(i;-,t). However, In the case 
where the plate is finite, the r formulation also contains z along 
the exiting plane GE (see Fig. 21(a)). As a result, the r and z 
formulations must be solved simultaneously. The surface tension model, 
described in Figure 22, also allows the solution of the r equation 
Independent of the z equations since z appears nowhere in the 
formulation. 

In any case, when solving the r equation, the finite difference 
representation of the problem results in N nonlinear algebraic equa- 
tions In N unknowns. A variety of methods were applied in order to 
obtain a solution to the simultaneous nonlinear equations. These In- 
cluded Lleberstein's extension of Youngs* work on over-relaxation to 
nonlinear elliptic partial differential equations (26), and the fa- 
mlllar Newton-Raphson method. None of the above methods were success- 
ful in obtaining a convergent solution. The technique suggested by 
Powell (3A) resulted In the method used in this paper to obtain solu- 
tions. Basically, Powell developed a subroutine wliieli was essentially 
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a "corapromloe between the Newtou-Raphaon algorithm and the methods of 
Steepest descents." In his paper, a Fortran subroutine Is described 
for solving the nonlinear set of equations, 

x^) 0 K - 1,2 N (138) 

Ttie objective Is to minimise the function 


K=1 


As with many iteration schemes, initial guesses are required for XI. 
This particular algorithm has an advantage in that the initial guessed 
values do not have to be that close to the exact solution. The com- 
puter program for the r solution contains the main program and three 
subroutines. The subroutine EQNS is the one supplied by Powell, The 
user supplies the subroutine MATINV which inverts the matrices and the 
subroutine CALFUN which contains the nonlinear functions F(Xi). A 
knowledge of the r solution makes the z formulation explicit in 
the unknown at each nodal point. As a result, a Gauss-Sledel 

linear iteration scheme was employed to obtain the solution. 

For details of the subroutines, the reader is referred to ref- 
erence 34. As we get into the computer results in the next section, 
a complete printout of the subroutine will be presented. 
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H Is chos.Mi to largo, l.e., „ choaon too large cr D ehoaen too 
small, all that is lost fs accuracy due to larger mesh npacings. For 
the conise mesh Infinite flat plate problem, the results of the numer- 
ical solution as well as the computer listing and final output can be 
found in Appendix D - "Computer Solutions/Listings." The solution to 
the r formulation (Fig. 20(a)) required 1098 calls of the subroutine 
used to solve the simultaneous nonlinear equation. The sum of the 

squares of the error to the exact solution was reduced to 0.00357 at 
the last iteration. 

The final values for the coarse mesh solution (r.z) were used to 
make Initial guesses for the values for the fine mesh solution. The r 
solution required the simultaneous solution of 159 nonlinear algebraic 
equations. The z equation, again being explicit in z^, resulted In 
164 unknown values of z^. The plot of the results from the computer 
program (details shown in Appendix D) are presented in Figure 24 in a 
print plot. The computer connected the nodal points with straight line 
segments. In general, there Is very little difference between the fine 
mesh and coarse mesh solution. The fine mesh solution required the ex- 
tended storage space option on the 1106 machine. The sura of the 
squares for the final r solution was reduced to 0.024 after approxi- 
mately 1000 calls of the subroutine CALFUN used to solve the slmulta- 
necs nonlinear equations. For this case 1000 iterations were required 
to obtain a satisfactory z solution. 

2. F. i iL i . te plate . - As mentioned in Section VII, Discussions of 

the finite plate formulation also contains z 
on the exiting jet surface (see Fig. 21(a)). 

of solution consisted of the following steps: 


As a result, Lhe method 
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Initially, aaauBad vnlun, of t along GE worn ohoaan. Iheno 
vnluoa worn naod to compute an t solution. The computed r solution 
used 161 calls of the subroutine CAhFim to reduce the sum of the 
squares of the residuals at the nodal points to 0.00066. With this r 
solution, the linear Iteration technique was used to calculate the com- 
plete 2 solution. Since the computed 2 solution resulted In re- 
fined approximations to the values of 2 along OE. changes could then 
be reflected in a new r computation. 

For the finite disk case In which the ratio of the radius of the 
liquid Jet to Che radius of the disk was one-half, a second r calcu- 
lation did not change when the 2 values were updated. A curve was 
faired through the calculated nodal points and is shown In Figure 25. 
Only two of the four available streamlines are shown In the figure. Wo 
attempt was made to refine the solution by completely doubling the num- 
ber Of vertical and horltontal grid lines. For this particular case, 
since no comparison with any existing analytical techniques existed, a 
more convenient value of 3. 3 was chosen for K*. The value of D em- 
ployed in the solution was 0.0138. The only specification along GE 
was that 2 was set equal to 0.25. The computer listing as well as 

the calculated r and 2 values at each nodal point can be found in 
Appendix D. 

The same method was used to numerically compute the finite plate 
case in which the ratio of the jet diameter to the disk diameter was 
three-fourths. There were 38 nodal points required for the r solu- 
tion and 45 for the z solution. A complete print-plot of the results 
Is shown in Figure 26. Two iterations were required for the z solu- 
tion as well as for the r solution. The 


sum of the squares for the 


T 
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r solution at the final Iteration was 0.001. The computer listing 
and printout can be found In Appendix D. Again, the value of K* was 
arbitrarily chosen as 3.3. 

The z solution corresponding to the r solution was obtained 
and is also presented In Appendix D along with its complete computer 
listing. The method applied to obtain the solution was a simple lin- 
ear iteration technique since equations explicit in can be de- 

rived both in the interior and along the boundary points. There were 
249 iterations required for the Z solution. Changes between the 
248th and 249th Iteration occurred in the fifth decimal point. 

A physical plane description of the infinite flat plate solution 
is presented in Figure 23. Only two of the four Internal streamlines 
are shown in the flgur-e. Curves were faired through the available 
calculated nodal points (r,z) by a best fit process. There appeared, 
at the onset, the question of whether or not the coarse mesh employed 
could sufficiently describe the flow. When applying boundary condi- 
tions at the free surface, a larger number of nodal points are de- 
sirable. As a result, the existing mesh was doubled. The total i|; 
was divided up into ten equal parts such that a = Ai|; = 0.05, and 
the total (|» was divided up into 16 equal parts. 

B. Surface Tension Dominated Model 

As previously mentioned, the r formulation for the surface ten- 
sion dominated flow can be solved independent of the z formulation 
(refer to Figs. 22(a) and (b)). The initial case examined had a Weber 
number of 4 and the ratio of the radius of the Jet to the radius of 
the disk was one-half. Also, the first solution to this problem as- 
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sumcd only axisymmetry . The results of the numerical soLutloc , epiploy- 
Ing a coarse mesh, resulted in a symmetrical r solution; symmetrical 
about the equipotential line emmanating from point D. This allnved 
us to make a simplification to the problem in that not only could axi~ 
symmetry be assumed, but also mirror image symmetry (i.e. , symmetry 
about the z 0 plane). This is significant for problems in which 
the surface tension effects are to be taken into account since the sur- 
face tension forces are highly dependent on the curvature of the free 
surface. By taking advantage of the symmetry Involved, additional 
nodal points can be placed on the free surface without using extended 
computer storage. 

Referring now to Figure 22(a), a vertical line was drawn (equi- 
potential line) emmanating from point D, where r - RD. RD is the 
dimensionless disk radius. The intersection of this equipotential line 
with the free surfaci^ AG was defined as point M. Along DM, it is 
known that z “ 0. In addition, the variation of r with i|; can be 
computed. Let us now refer to Figure 27 in which the equipotential 
line DM is depicted. Since DM is an equipotential line, the ve- 
locity along this line must be equal to V, the incoming Jet velocity 
with the exception of the point r = L in the physical plane. At 
point D, a velocity discontinuity will exist. However, we can specify 
how r varies with along this line as follows; 

In the physical plane, 

On DM 

u=0 on z=0, Lir^R (140) 

’ max 

V == -V 

where R is the maximum radius of the liquid flow pattern, 
ma:" ^ 


J 


L 
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An expression for the stream function along DM can now be de- 
rived since the radial velocity component along DM is 0, we have 

0 - - -i (1A3 

r 3z 

This implies that t|» = i|>(r) alone. In order to find what the function 
is, the definition of the axial velocity is employed, namely, 

V - (l/r)(3iir/3r), on DM 

-V - i (u: 

r 3r 

Integrating this 


+ Constant 


Applying the boundary condition that iJj = 0 at r = L, the constant in 
(143) can be calculated. The following expression results, 

^ . Ve! + Vk! (144) 

If this equation is nondlmenslonallzed. 


■ 2 ' 


As we have done in the pist derivations, the starred notation is 


dropped. 


- - (-] 
2 2 J 


Solving for r, 


r = +\ 2 \\> + (RD)‘ 


where 
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ECaa.lon (147) represents the bounder, condition en^loyed nlon* DM In 
the inverne plane tor the r formulation. The condition t - 0 wan 
usi'd for the z formulation. 

In the courne of finding the solution to the r formulation, de- 
picted in Figure 22(b). It was necessary to solve a cubic equation for 
the parameter T along the free surface. Physical as well as mthe- 
matlcal interpretation was required when choosing the proper root of 
the cubic since a poselblllty of three real root, existed. Referring 
to equation (C-78) In Appendix C. the only mathematically meaningful 
roots are those In which the absolute value of the parameter T was 
greater than or equal to one. However, the possibility still existed 
that all three roots would be real and In addition satisfy the re- 
quirement that their absolute values were greater than one. As a re- 


sult. some physical insight was required when deciding upon which roots 
to employ in the equation relating T to the fictitious point f. 

(eq. (C-84)). For example, it Is known that as the nozzle exit is ap- 
proached along the free surface, z = f(r) becomes steeper (l.e.. 
f(r) approaches Infinity). This would coincide with the curvature 
terms dropping out of the boundary condition In the physical plane 
formulation. An alternate approach to viewing this is that the fic- 
titious point t approaches the Image point z^. As a result, the 

value of the parameter T approaches unity. The algorithm selected 

for choosing the proper root of the cubic was to select the value of T 
closest to unity but ensuring that Its absolute value was greater than 
or equal to one. During the course of finding the solution, problems 
in implementing this algorithm occurred. parUcularly when close to 

nozzle, since the T values closest to one wore slightly less than 
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one and wore automatically discarded by the algorithm. The resulting 
potential lines and atieaialinos appeared inaccurate when plotted up. 

The only way found to circumvent this problem was to set up an addi- 
tional algorithm which set T identically equal to one for several 
free surface nodal points in the vicinity of the nozzle (l.e., those 
nodal points in which the computed r value was <1.00 N). Physically, 
this reasoning is justified since it is known that f'(r) must approach 
infinity there. 

The results of the numerical solution are Indicated in Figure 28. 
The Weber number for this solution was 4 and the ratio of the radius of 
the liquid jet to the radius of the disk was 1/2. A value of 3.5 was 
chosen for K(j>, and D was set equal to 0.0204. The computer listing 
can be found in Appendix D along with the computed r and z values 

for the nodal points. The sum of the squares for the r solution was 

“6 

0.77x10 and since the r solution was independent of the z solu- 
tion, a second iteration was unnecessary. 

C. Discussion of Zero Gravity Results 

The numerical solutions were compared with the available semi- 
analytical results of Schach (39) for the case of the infinite flat 
plate. The method employed by Schach is attributed to Trefftz. The 
method used by Schach Jld not appear readily extendable to more com- 
plicated geometrical flows and could not be employed to account for 
the effects of surface tension. In making the comparison between rel - 
erence 39 and the numerical results, the fine mesh solution presented 
In Figure 24 was used. The results of the comparison for tJie liifJnitc 
flat plate are shown In Figure 29. The symbols Indicnled in the figure 
were obtained from Figure 11 of Sclinch’s paper by using an exnaiulnli 1 <■ 



f)4 

„calo. A, a raault .hara ia »o„a ualcna^ arror aa.oclaaad wie, tUa 
procGSa of taking the vesulta from the reference. In any case, 

.,raa«,„t looR. partlaularl, «eU wUh the sola excaptlaa of U.a firac 
. coordlnata greater than unUy. C„a fioal podat with raopact to tha 
Infinite fUt plate eolutlon conoetna the extretne left coordinate 
Figure 29. naaely. r - 4. a - 0.125. Thane two values are fixed hy 
continuity, both in our numerical program and in the aeni-analytical 
neaulta of Schach. Thin result was ohtained as follows, a..u.ing con- 
ntant density, the volumetric flow rate into the control volume at A3 

coordinates, the flow is given as and the flow out hy 

2nR. ZnV. R,uatlng these and cancelling leads to the fact that 2„ 

musrloual R^/dR^,,)- Nondimensionaliting with respect to R„ yields 

i (1A9) 


* » - 

■ 2R* 


jet 

Dropping the starred notation and obseiv/lns that R^^^ ' 4 at the 1 
houndary of the control volume shows that 2„ must equal 1/8. 

AS far as the finite plate is concerned, there was no avallahle 

comparisons with past experimental or analytical worh. As a result, 
comparisons were made with respect to our own xero gravity expcrl- 
aental data. The results are shown in dimensionless coordinates in 
Figures 30 and 31. Figure 30 is for the case whore the ratio of the 
Jnc radius to the dish radius is one-half and Figure 3. Indicates the 
comparison when the ratio of the let radius to the disk radius is 
H, co-fourths. The comparisons were made with respect to the outer or 
„,p free surface since it was impossihie to view thf >ower free surface 
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because of the way in which the flow occurred. The results were gen- 
erally good for both ratios compared. Experimental data points were 
obtained from both sides of the axisymmetric sheet as it flowed around 
the disk. An averaging procedure was subsequently used to plot the 
continuous lines indicated in Figures 30 and 31. The analysis corro- 
borated the experiments in the sense that as the ratio (Rq/I*) becomes 
siviller, the jet appears to leave the disk more tangentially. 


IX. NORMAL GRAVITY EXPERIMENT SECTION 
A. Apparatus and Procedure 

1. - The experimental teat rlj need te obtain the 

normal gravity data la shown In Figure. 32(a) and (h). The rig con- 
elated 01 an angle-iron frame In which wan mounted a 10 gallon euppl, 
tanh, a settling chamber, a 56 gallon catch basin, a return pump, a 
control box. a cloth, aeguence timers, a regulator, and a supply tanh. 
The major lunctlons were controlled through the control box. * high- 
speed Mitchell Monitor motion picture camera (nominal speed 
400 irames/sec, was located directly In front of the experimental 
test rig. The camera was mounted on a IPsllensah camera stand which 
was fastened to a concrete floor by means of conduit clamps. 

The experiment could be conducted In either a pressurised or non- 
presaurlsed mode (gravity-feed). To operate In the pressurised mode, 
two vent valves located above the supply tanh were closed and the sys- 
tem was pressurised through the regulator. The pressure level was re- 
corded on the gag. located l^sedlately to the right of the regulator 
as shown in Figure 32(a). For both modes, the return pump was used to 
resupply llguld to the supply tanh In order to maintain a nearly con- 
stant level of llguld. Since the supply tanh was fabricated from 
stainless steel and provided no visible means for monitoring llguld 
level, attempt, were ,»de to connect a plastic hose between the needle 
velve located Just upstream of the solenoid valve and the side of the 
supply tanh. This system was generally Inaccurate and as a rule the 
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screen was mounted below the f^r=^ 

the first screen end both screens were butted 
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Of the spacer. This tschnlque ssened to elim- 
inate sn, noticeable swirl 1 „ the flow. 
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etliniiol and dlutlHed water. Tholr proportlea at 20“ c are listed In 

Table 1. Nn attempt wan made to eortect tho.,e fluid properties tor 
tempcraturfi changes, 

1. Tjat.j.n,eedu«. - Prior to „ test run the settling ehamher was 
fUltd with tl.e teat liquid by opening the solenoid value while holding 
the bottom of the settling chamber closed. The liquid completely 
filled the settling chamber until It flowed out of the relief screw. 

At that point, the relief screw was closed and the outlet could be 
opened without loss of liquid from the settling chamber since it was a 
stable pressure supported system. This method worked when using dis- 
tilled water for all the nestle openings. It did not work for anhy- 
drous ethanol in conjunction with the largest or 1 centimeter diameter 
opening. As a result, for those tests, the nettle opening was kept 
sealed until the solenoid was opened at the start of each test. 

Calibration tests were made prior to every series of teat runs. 
Electrical sequence times In Che control box were used along with 
Ltadcated cylinders to determine the volumetric flow rate. At least 
three calibration tests were made before each series of runs, and an 
average value for the flow rate was thus determined. The Impingement 
velocity (velocity at the disk) was calculated by correcting for the 
effect of gravitational forces. 

8. Results 

1. SJM!lyrstatR_nowjm - Basically, two types of flow 

patterns were observed 1„ the course of normal gravity liquid Jet 
Impingement, The first type occurred when „ very high speed Jet im- 
pinged upon a nelld aurfaee and spread nut rndinlly as a thin aheet. 

“"‘'1 it l»>t..me naslable and broke np late 


L 




Rmall liquid droplcits. A second pntU'rn Is shown sc.homall rally In Flp,'- 

uio :)5. Ko bioakup was obs. rvod In t.lilr; pattern; tho jot: rnrvod upon 

lusc'lf to fonn a surface of revolution or plunit*i The Incoming, .jot vo'- 

locity Is V and tlu,' radius of the circular liquid jet Is The 

disk radius Is denoted ar, L whllt^ H Is the distantly between Chf* 

no?.zle and the disk, Tlte maximum radius that the plume prsse.ssed *s 

denoted as R . The maximum plume radius war. the prluviry experiineiU.;] 

P 

variable for the normal gravity rtudy since if was easily measurable 
and characteristic of the entire flow pattern. For the case when the 
plume was not visibly disturbed* the plume would slowly move to a 

final R , with visible surface ripples. 

P 

2. Unstable lumps . - In attempting to obtain steady-state Im- 
pingement profiles, an unusual phenomenon was discovered. Instead of 
a single steady-state flow pattern, a number ot unstable flow patterns 
were observed prior to the attaining of a stable steady state. Typ- 
ically, the jet impinged upon the solid and formed a plume of some 
given radius with surface ripples present. If disturbed, either up- 
stream or downstream of the plume, the jet would jump to another ap- 
parently stable configuration of larger plume radius. Normally, only 
one jump would occur in a single test run, but occasionally rwo jumps 
or three apparently stable configurations were observed. A series of 
teats indicating the flow patterns before .and after a jump are shown 
in Figure 36. Figures 36(a), (c), and (e) illustr.ate the beh.avior be- 
fore the respective lumps, while Figures 36(b), (d) , and (f) occur 
after the jumps, The plume size for the final conf Igur.itlons appear 
to 1-: nearly double what they were Initially, However, there does not 
appear to be any correlation between Initial and final size, Tlie jumps 
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«ft U,c ,H.,k „„,,|,r ,,„„ ,„„„. „„„ ,,, 

inltlntod by ,. dj «,.„rbm,n,.. „„ „ u.t. fj„„, 

TIu, «c,„,|y-„f„t„ «„.flgurbtln„ J, ,;„ 

turbances, l.a., u at aomu tlxad plum,- aba,.,.. u„d„r „„ ,.|r- 

cimatancea ™a tha pluae obaoryad to Ju»p u, a am,-, II,., plom,. tonrig. 
oration: 'ill .lumps vara to larger plume altoH. 

The actual causa of cheae Jumps remains ima„soar,,d at tbl.s writ- 
ing. Soma possible aauaea have been eliminated, such as awlrling fl„„. 
t£ the flow were swirling, the ateady-atate flow patterns would be 
aigniflcntly affected if half the flow over the dlak was obstructed. 
However, the pattern was not affected at all when this was done. It 
la tentatively concluded that the ,1umpo are natural in occurrence. 

The initial states and all transition states are unstable to small 
disturbances. For the purpose of the experimental Investigation, then, 
the steady-state was chosen as the state in which futthe, dlstnrhaaccs 
caused no change in the liquid flow pattern, 

3* KxEe r i jnental data . ~ The experimental runs were conducted .such 
that the viscous dependence would be small. As shown 1„ section HI, 
depending upon the ratio . there exists a crltloai Reynolds nnl- 

hor above which the flow can b. consldared as essentially Invlscld. 

The Reynolds number was calculated at the point of Impingement on the 
disk (i.e., the approeoh velocity to the disk was corrected for gr.iv- 
Itatlonal effects). From analytical oonslderatlona, numerous dimen- 
sionless parameters arise and thus form the basis for cortel.iling , |„. 
primary variable, the plume site. These parameters l„cl„,le th„ „„h„.. 
nunber, which Is the ratio of inertial to surface .ension fnrees. ,u„l 
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the Bond munber, which is the ratio between nravitatlonnl and Htirfare 
tension forces, In addition, several geometrical ratios appear, such 
as Rjj/L, the ratio of Jet to the disk radius, and H/R^, the ratio of 
the nozzle height to the jet radius. The experimental results arc 
listed In tabular form in 'fable 3, The last column contains the plume 
radius, nondimensionalized with respect to the jet radius. In ex- 
amining the difference between zero and normal gravity liquid jet Im- 
plngemert, it can be seen that two additional parameters are required 
in order to completely define this phenomena in normal gravity. They 
are the Bond number, Bo, and the dimensionless nozzle height, H/R^. 

The Bond number relates the relative contribution of gravitational 
forces, while H/R^ has an effect in that jets flowing downward under 
the effect of gravity accelerate and thus shrink in size, thereby hav- 
ing an effect on the resulting flow. 

Data analysi s. - A linear regression analysis was employed in 
order to correlate the independent variable, Rp/Rg, with the remaining 
system parameters. Let us define the following variables, 


R 



(150) 

0 

B 

Z = 

1 L 

(151) 

z 

^2 R 

(152) 

o 

= We 

(153) 

Z, = Re 

(15/.) 

Zj = Bo 

(155) 

The parameter 

Z^, the Reynolds number, was inclutU'd to see wli.H 
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the viscous Influence was Th« 
following complex model 
R 

^ = 22.8241 + 7.67190 Y „ 9 

o /• 0/129 X, - 41.2001 + 42 , 

^ ^1 ^ "»2.9330 Xj^ + 1,91085 X^ 

+ 8.90368 X 3 9 . 1.46319 X, - 2.07823 X, * 0.0696939 x.x 


(151 


~if ' 0.371839 
^154239 


H 


^ ~ 9.22727 


2 Tag^i— ' 

X = We - 129.227 

3 75.7108 

X, = Ro - 4663 .ro 
^ 1897.57 

X^ = Bo - 2.15602 
0 1.12058 


The value of r 2 above statlsf-f 

auuve statistical model wa<? n qbi 

"•sans that 98.1 percent nf tn 

cant of the total variance In g /n .. 

for bx the model An .. ° ^ accounted 

model. An examination of the al.„e i„ ... 

complex model telle us ho» R /p . 

R /R f ° ™ ® “f fha parameters- 

p'®o Increases slightly as h/r end a 

^ o are increased, and in 

creases significantly with x tb u . 

vo.vcd in t„o tetm, l„ th 

-angest tetm (1 .. -apenda on the 

suit. '-“"''""’^'"'-^--fflcleno. Asetc- 

»p/R„ dectMae,. with Increasing Bo n,e „ . . 

"oru canr„.,i„g sine.., u „ , . ' *>= 

“ -rms. The actus, data 
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bears this out, R /R sometimes Increasing, or sometimes decreasing, 
p o 

The coefficient for la essentially the largest of all (X^ Is re- 

lated to We), kfj a result, a simpler model was pursued using only Wo 
as the Independent variable. The following equation resulted; 

^ « 21.7059 4- 10.1410 (lfi2) 

o 

2 

For this simple model R = 0.916, which means that the correlation 
accounts for 91.6 percent of the variance in the data, A plot of the 
data Is displayed in Figure 37. 

The Weber number turns out to be the most statistically signifi- 
cant variable. This Is not to say that the other variables are not 
significant, but that the Weber number Is the most dominant variable 
in determining R /R . 


X. CONCLUSIONS 

- An nxporl»nntal and annlyUnnl Inva.,tlg„,l„„ 
nondunted t» dntorn,l„n the free surface shape of circular jcta 

Impinging uuraal to sharp-edged disks 1 „ aero gravity. The test 11 ,- 
ulda etiployed were anhydrous ethanol and trlchlorotrlfluoroethane. 
dot radii uere varied fro. 0.15 to 0.75 centloeter and disk radii of 
1.0 and 1.5 centlaeters were eaployed. The Jet velocity was varied 

n 12 and 365 centimeters per second. Under the stipulation that 

the nozzle was located lonoi- c . 

least 5 centimeters from the disk, the in- 
vestigation yielded the following results: 

1. It was analytically deteralned that there exist flow regions 

fhere viscous forces are not significant when computing free surface 

*Napes. It was shown that the Reynolds number pVR„/u and the ratio 

»f Jet to disk radius R„/L uniquely define the flow regions. It 

»es further shown that the Reynolds number specifying the transition 

between viscous and nonvlscous flow decreased with Increasing Jet to 
disk radius ratio. 

2. Within the invlscld region, three distinct flow regimes were 
experimentally found which depend uniquely on the kebor number 
. V R„/o and the ratio of the Jet to disk radius R„/,„ These flows 
eflned as Surface Pension Plow, Transition Flow, .and Inerr in 
Mow. The critic.,, „,.ber number between rcglm..s was found to 
with Inrrca.slng U-t to radiu.s nuio. 

'■nio surrof... 
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lino. »a. obtainod for tho ca.o of imploRomon, noro.nl to an Intlnltn 
flat Plato and comparod favorably with .aoml-analytical to<lmf<|oon In 

the literature. 

4. A naoorlcal solution yioldlng froo aurfaco shape, and streao- 
llno. wa. obtained for Inortlally dominated flows at ratios of Jet to 
disk radio, of one-bait and threo-fourths. The comparison with ex- 
periments showed good agreement for the upper froo surface. 

5. A surface tension dominated flow was formulated and solv, , 
numerically. The system Weber number was 4.0 and the ratio of the let 

to the disk radius was one-half. 

normal eravlty . - An experimental Investigation waa conducted to 
determine the characterietles of circular llguld Jets Impinging nor- 
mal to a shatp-adged disk In normal gravity. The te.st liquids em- 
ployed were distilled water and anhydrous ethanol, dot radii hetweea 

0.125 and 0.50 cantlmater wera employad and the disk radii were varied 
between 0.5 and 2.0 centimeters. The Jet velocity had the range of 
75.5 to 484 centlmatera per second. The distance between the ncsxlc 
and dlak was varied between 0.25 and 5.0 centimeters. Under the 
atlpuUtlofi that the Reynolds numbers were such that they exceeded tin 
minima value required to avoid viscous Influence, the investigation 

yielded the following results: 

1. The liquid flow pattern was observed to Jump from uue .,p- 
parently stable flow pattern to another until .1 comp,e.e,v ....ihic c,„.- 

flgnratlon was reached. The Jumps wore triggered hy i..u h 

upstream and duwn.troam of the disk .and were uppnrentl, u.U,„.al - 


currencG. 


2. The dlmenslonles. plume radius R/R, w,„. , err. , u. e.l by 


i 


76 


of a linear regreusion analysis . A slraplc model, employing only the 
Weber number, accounted for nearly 92 percent oi the experimental data 
The following empirical formula resulted 

% We - 129.227 

^ - 21.7059 + 10.140 —5.7108 

0 

where R is the plume radius, R the nozzle radius, and We the 
P ° 

system Weber number. 


XI. API'KNDIXF.S 


Appendix A -Detailed Cnlculationa of Free Surfneo Boundary Condi Lions 

There are two boundary conditions required for tl>e case of a free- 
surface in a fluid dynamics problem. This is In comparison to known 
boundaries in which only one boundary condition is required. Tlie two 
conditions to be ratisfled are: 

(1) Conservation of energy along a streamline 

(2) The velocity normal to the streamline is zero. 

Consider the following geometry: 



Figure A. 1. - Schematic of Liquid Jet Impingement 


where "i" represents any point on the free surface = f (r) and "d" 
designates the reference point which Is chosen as the point where ttu- 
liquid jet exits from the nozzle. 

(1) Conservation of energy along n stTnanllnc* 

Bernoulli's equation .rlttcn between points "i" and "d" be. 


1 

2 


2 2 
( II + V ■ ) H- 




P, 

1 


on V 


) ( t- ) 


(AH 


8 


The pressure at point 1, p^, Is not the same as it Is at point d, 
Pj, due to the effects of surface tension. In general, 

Pg - P - JJ (A-2) 

where Pg ° known gas pressure, and 


«2 


where radii of curvature where 


r 




and 


(A- 3) 


(A-4) 


^2 / 2 \^^^ 
r(l + f'^j 

Combining (A-^) and (A-5) with (A-2) we find 


Pg - P « 0 


f" 


f 


L(x . vf' 4 . rf 


which can be combined to yield, 
Pg - P 


0 d / rf \ 

Vi + f'v 


Applying at point "d" 

2 = H, r - R^, and 
For large f * , 

}jl + f*^ A ^7^ = f» 


df . f, 
dr 


(A-5) 


(A-6) 


(A-7) 


Therefore, substitution into (A-7) yields 


Pg - Pd = 


0 


(A-81 


/9 


ApplyliiR at point 


Pg - Pi - 


rf ' 


^ ^ 1^1 + f ’7 

Now (A-8) and (A-9) con be aubstltuted into (A-1) to obtain 

' rf \ 


1 / 2 . 2v 0 d 
Y <u^ + V ) - — 


i v2 - 


pr dr ^ ^ f,2j 2 ■ pR 


on Ky - f(r) 


(2) The velocity normal to the free surface Is zero 

A portion of the free surface z = f(r) is shown in Figure A 

s 



Figure A. 2. ** Velocity Vector at Free Surface 

n is the unit outward normal to the surface at some point and u 
V are the. velocity components such that V = ui + vj . Since the 
loclty normal to the surface is zero 
V • n = 0 

With the surface given by z = f(r), the unit normal is given as 

s 


yi + f '2 


Hence 


> - uf 

V • n 




yr+f '2 yi + f 2 

which simplifies to 


= 0 




(A-9) 

(A-10) 

.2 

and 

ve- 

(A-11) 

(A-12) 

(A-1 3) 


-uf + V = 0 


on 


(A-J4) 


Appendix B - Zero Gravity Drop Tower Test Facility 

The experimental data for this study were obtained In the Lewis 
Research Center's 2.2~Sccond Zero Gravity Facility. A schematic dia- 
gram of this facility Is shown In Figure B.l. The facility consists 
of a building 6.4 meters square by 30.5 meters tall. Contained within 
the building is a drop area 27 meters long with a cross section 1.5 by 
2.75 meters. 

The service building has a shop and service area, a calibration 
room, and a controlled environment room. Those components of the ex- 
periment that required special handling were prepared in the con- 
trolled environment room of the facility. This air-conditioned and 
filtered room (shown in Fig. B.2) contains an ultrasonic clv^anlng 
system and the laboratory equipment necessary for handling test Mq- 
ulds. 

Mode of operation - A 2.2-second period of weightlessness is ob- 
tained by allowing the experiment package to free fall from the top 
of the drop area. In order to minimize drag on the experiment package, 
it is enclosed in a drag shield designed with a high ratio of weight to 
frontal area and a lov drag coefficient. The relative motion of the 
experiment package with respect to the drag shield during a test is 
shown in Figure B.3. Throughout the test, the experiment package and 
drag shield fall freely and Independently of each other; that Is, no 
guide wires, electrical lines, etc., are connected to either. There- 
fore, the only force acting on the freely falling experiment package 
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r 


I 


Which is estimated to be be3ow lO'^ «•«, 

Mease,^y^. . ^he experlmeat package, Installed within the 
ra, ahlela. ,, ^ 

*hX. atraaae. ,„p,e „,,p 

aaaa syatam ooasista of » do„bla-ooti„g air cylinder with a hard- 
atael hnlfe ed,e attached to the piston. Preaaurleatlon of the air 
cy x..aer drives the hnlfe edde against the wire, which la hached hy 
a. anvil. The reanltlnd not.h caneea the wire to fall, snoothly re- 

aaalnd the eaperlnent. No n^serahle dlatnrhances are leperted to 
the package by this releeee procedure. 

Mcoverv svat.». . gfr„ ^he erperlment package and drag shield 
-e tr™ the total length of the drop area, they are recovered 
ace oration In a 2.2-e.eter-deep container filled with sand. The 
-aaleratlon rate leveraging IS g-s, is controlled by selectively vary- 
ing the tips Of the deceleration splkea bunted on the botton of the 

g shield (Pig. B.l). At the tine of Impact of the drag shield In 

the decelerator container, the exDer^^,d.n^ 

*P ent package has traversed the 

vertical distance within the drag t*hielH ( 

arag shield (coc[.are Figs. B3(a) 

and (c))e 


Apiu'luilx f' - 0('t a / U-<| Dnriv.iUniiH of Kin I to l)iff<Toiu'o 
OporatOTH Alonp JUminlar Ion 
A. Infill I to Klat Plato 

Tho flnltu cllfforonoo roproanitat Ions for tlio dorivat I von ai 
llioao shown in the toxt (oqs. (U)) to (liO)). kofor to Klpiiron 
I '3(a) and (h). 




(C-l) 


Equation (C-1) bticonic'H 


«■) 


r - r r„ T r -v 
o o 2 o 2 
a 


Or 

fiip 


0 on CE 


-~i "'i '■ "i* ♦ A ‘■'i - '^y 

/ • Off 


V 



Application of 3r/3ijj on CE yields the fact tiiat r., " r, 

fc. H 

does not Involve the unknown fictitious point f. However, 
Implies r ^ Therefore, r must only be a function of 

(j>, r = But, along GE, (t = constant, (iji =• 0). Hence, 

r = constant along GE. 




The r difference equation can be written 


I 


(•: <i) 


and 
= 0 


»!/» 


-I - I <>f 


] 1 


(r 
H I 




2f( 




. (r, - r,,,' 

K(» 


(('.-'>) 


MomR AC‘- 


(to 

This Iniplips 


2 w„_v2 


+ r 


2/3r 


\:j4 + m' 

2 y o\ 2 ) 


Rearranging and letting Aiji 


a (recalling a ~ Av/ Av) 


(i-f - 2 


2 ^ 1 
r 

o 


4=2 - i C) - '^3>' 


\fli ere upon we can 


calculate the two expressions 


r, - r, « - '3>1 


Tf + 


2 ^’'l " ’’3^ 


+ 2r, 


(oi.) 


(OV) 


(C-8) =nd (C-« into (C-5) ylnldn tl'o rvlntl™ 

1 ''o 


(C-H) 


(r-^0 


^ A, . rf, |i - i - j. (t, •■ n/] r 


2, (r, - 4,. 


Lu 


(r, -V r.,) 


.1 0 


(c- ln> 


'/ « Constant on AH, let x " " 


j 
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The z difference equation can be written 



which can be expressed as 



This yields the two relations. 



(C-17) 


(C-18) 


(C-19) 


(C-20) 



Inserting these last two equations in (C-17) eliminates the fictitious 
point Zf, yielding 




o 






r a 
o 


Applying the boundary condition along 


+ 1 

AG (see eq. 


(C-20)) 



but ^2 “ 

h- 

Equation 



z^, hence we write 



(C-23) becomes 




( 025 ) 


(C-26) 


^ ^ ’'oJ ^1 (C'27 

’‘o ■ 2 1 + 1 

+ “ o 

such that point G. 

B. Finite Plate 

Referring to Figures 16(e) eud (b). the chengee in the boundary 
conditions between the Inllnlte end finite plate occur on GE end the 
addition of the free surface EG. In edditlon. G become. . speclel 
point in the r formuletlon while both G and E become speclel 
points in the z formulation. 

- 0 on GE 
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L 


’© 


i© 


Along GE, the r difference 
_3 


equation la 


r - 


T 


• v • '6 - s '■= - V'] ■ A ... . 


f) 




in finite difference for., the boundary 


- f r 


condition la, 


A'* Aifi ^^2 ~ ^4^ 0 

solving for f, 

i “ r + ctr (z^ - ss ^ 

^ o' 2 ^4^ 


(C-28) 


(C-29) 


(C-30) 


and, inserting this Into the r dtff 
result, ° tho desired 


''o ■ f <<:■> + r, ) + - i . ’1 

oU ■ « - -2-i^^i - 




* 1 


oil ED 


o / 2 

a (^2 “ ^4)^ “ 0 


(CN31) 


I 





The r difference equation can be written 


-O - T '^2 ^ '£> 'o: 




+ r^) 


+ "i - ' 3 ’^ " “ 

8a 


(II) + =1 


+ r^ == 1 

i 2 AA I o\ 2 AiJ) / 


Yields 


(r - r )2 . ip - -\ (^ - '3>^ 
r L a J 


and 


(C-32) 


(C-33) 


(C-34) 


r2 + Tf - 2r2 - \ ~ ’'3^ ] 

Substituting (C-34) and (C-35) into (C-32) we obtain 

- ^>2 ^ 'cfi " i VP ■ ' ''fl ■ ' 

- * \ <'•1 - ' " 


(C-35) 


r. 3 ) 


(r-30 


1 


Now, points K and 




puintH Minn.. nquatJon (f~H) 



Applying equation (C-10) at point G, 


'o - 'o-4 + -I {^2 - |^4a^ - A (r, + ,p 


2 a 

® f? = 0 


(C-37) 


Now, applying + r/'-^V^U 


HJ \3t/i 


fh - ^ 2 

1 2 A.f 


^j.j- 0 along GE without Involving 

(C-38) 


- r,\/z. - z 


+ r ( --- ^ \( ° ~ 4 

“I A^'“” 


Solving for f^, 

fj = rj + 2„r^(3^ - 

Therefore, equation (C-37) becomes 


(C-39) 


■*' [*'o^''o ■ "'4^^ ■ tJ ” ° (C 


V^V'/ U^Adr/ “ ^ "" 


O'j 


Application of z difference equation yields, 


r "t* f (^2 ** 

J 1 l\ 2 fr^a^ + l) A(r^a + 

^ (C-«) 


Now, along GE, r^||^ - (if)(S) “ difference form can be 

written. 


'■■ - (=H^xp-;) ■ " 


'o\ 2 AiJ» 


Solving for f , 

f = - r^«(r2 - r^) 

Hence, equation (C-41) becomes 


r 2z, - r a(r. - r, ) 

o / . N . 1 o 2 4 


(^2 - 


/ 2 . 1 \ 

(r a + — 

V o «/ 




i 


J 


© . 


0 1 


The* V. <1 J f forcncc c(]uni:Jon rni) bc> wriU-tn, 




(z^ + Z^) 


/ 2 l\ ^^2 ~ 

nv^i) 


2 ? 

. T , 2/3z\ . /3z\ 

Applying r / — 1 + f— j = 1 on ED, 


^2 / ^ . 2 . - ’’■f f ^ / h ~ 1 

O I 2 A«i- I I 2 A(^ / 


Solving for (z^ ~ 


"2 ' "f 


= -I - i u, - . 3 )^] 


Also , 


+ -^2 " 2==2 - J - -L . ,2 

OIL n J 


(C-45) 


(C-46) 


(C-47) 


(C-48) 


Substitution of those last two equations into eq’jition (C-43) yields the 
clcsirod expression 


I 
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Applying equation (C“22) at point G yields. 




Solving for 

r„ - , - 2r u(r - r . ) 

2 s. n o 4 




(C-50) 


(C-51) 


( 052 ) 


Substitution into (G--50) yiolds 



In ordor to find the fictitious point, ‘*PP^y 



ca 


0 on GE, wltlvmit liivolvlnp, f 






lift! V lit}' for 1^, 


( 


7 




Finally, oltrain I ii>', I roin (‘qua t Jon (0')4), 



Applying r difference equation (eq. at point o when.' f 

is a fictitious point outside Ihe boundary 


f 


1 


3 

17 

o ?. i 4 o 


[„2 - J <-^t - - J 2 <'l 


+ r^) 


.^(rj-r/.O 

in difference fornix the free surface boundary condition becomes 


IL1I£\ m 2/5iJlIiV 

2 j ^o\ 2 I 


1 + 

We We r dr 
o 0 


r^)/2 Aif 2 
r^)/2 Aii» 


Simplifying by using a = A<|t/ArJ» and defining a “ AiJ; yields, 


-i (-^1 - -^ 3 )' * - v' - 

a 


1 2 2 1 d / ^0 

-we We r^ dr 


Now, examine 


_L -L. 

r dr 
o o 




As an approximation to this derivative Q is assumed as a constant. 
In actuality, Q *= f(r^) and r^ » Expanding 


1 d / ^o 
^o ‘*’^0 I yq 4- r2 


yields 


f 


I 


A -A 

r dr 
c 
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2Q + r; 




Therefore, the finite difference representation for the free surface be 


(C-61) 


comes 


■\ (r^ - r^)^ + r^(rj - r^) 



(C-62) 


Tj must be eliminated between equations (C-58) and (C-62) . Let us de- 


fine 

X - “ ’'a 

Equation (C-62) can be written 


(C-63) 


1 / ^2 . ^2 2 _ 
~z <’'1 ■ '' 3 ’ * V ■ 


4a‘ 


2 2 
^ We We 


2Q + r;: 


y<! + + 'o), 


Where 


(C-64) 


„ - -3>. ■ _L 

Inserting (C-63) into (C-58) yields 


(C-65'» 


2 2 ,3 

V 


2 

. r *■ 

4 L 3 o L o 

r + r-’r . r H 

o o 4 ^2 „ 

a 2a 




Tl'is is of the form 
'X 


A'x^ + B*x + C = 0 


Hence 
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-B* + Vb'2 - 4A'C* 
X 2A' 


(C-67) 


2 

r 

A' = — 
* 8 


(C-68) 


B' = -;r 


2 

C* •= -r^ + ? + — % 2 


(C-69) 


(C-70) 


The sign In front of the square root in equation (C-67) was chosen as 
positive since x greater than or equal to zero. In addition, 


>/3z\^ . /8z\ 
\ai|i/ \ h ) 


^ We We r dr 




on AG 


Applying z difference equation (132) at point 0 where f is a 


fictitious point outside the boundary 




if + z.) + 




r (zj^ + Z3) + 


(f - Z^)(Zj^ - 

a) ' 
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In finite difference form, the free surface boundary condition can be 
written. 


^ <'l - “s* 
a 


where, 



(C-72) 


Q* 


2 " V 

l " ■ " ' ; 


(C-73) 


As an approximation, Q* is assumed to be a constant. Expanding 

\ 



V^oQ* + V 1 


(C-74) 


Therefore, the finite difference representation along the free surface 
becomes 


r^ff - z.l^ + -4 fz. - z ^ 4a? ^ 

""o (r^Q* + l)yr?Q* + 1 

(C-75) 

f must now be eliminated between equation (C-71) and (C-75). The man- 
ner in which this is done is as follows; Equation (C-73) is solved for 
f ® ^ct^^l'^3’*4^ * results are then Inserted into equation (C-71). 

In this way, equation (C-71) remains explicit in z^. Actually, it 
will be more convenient to solve for the variable (f - z^) instead of 



10? 


f since* equation (C-71) can be expressed as, 


— , 2ar‘ + (z ~ z ) 


o 2 
o 


r2 + 4 


Turning our attention to equation (C-74) , it is solved for (f - z ) 


(-■l - 2,)^ 


+ 1 a 


4a^g^ 


1 - J- + ^. X 

] We We r f 2 2>^ ,2 

° vii~ 


(Zl - Zg)' 




A new variable ig introduced 


(C-77) 


'*1 - 


(C-78) 


Also, let 


<'1 ■ 


(C-79) 


C = 1 - ~ 
1 We 


(C-80) 


C. = — • 

2 We r 

o 


(C-81) 


Making these substitutions into equation (C-77) results In a cubic 


equation foi T 
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Rewriting this as follows; 


(C-82) 


+ aT + b = 0 

Once T is found from (C~83) , (f - z^) is calculated from equa~ 
tion (C-78) as follows: 



o 


(C-83) 


(C-84) 


A)>])crid Jii I) • "(lonipiit f’i- ,So1nl ioiis/Lf 


A. InliiiiLo Flat Plate 

1. Computer Ll.-Jtlng 

11. Coarse Mesli Solution for r and ?. 

111. Computer Listing 

IV. Fine Mesh Solution for r and z 

B. Finite Plate 

V. Computer Listing (Rq/L “ 1/2) 

VI. Coarse Mesh Solution 

VII. Computer Listing (R^/L = 3/4) 

VIII. Coarse Mesh Solution 

C. Surface Tension Model 

IX. Computer Listing 

X. r Solution 

XI. z Solution 
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Figure 6. - Flow system schematic. 
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Figure 10. * Zero gravity experimental results. 
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Figure 14. - Dimensionless governing equations and boundary conditions 
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Figure 16 (a). - Inverse formulation excluding surface tension (r solution) for finite plate. 
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Figure 16(f)) ■ Inverse formulation exclutling surface tension (/ solution) 
for finite plate. 


















Figure 19. - Nodal point representation for rec- 
tangular mesh. 
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Figure 20(a). - Finite difference representation for infinite piate (r solution). 
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gure 31. - Comparison of numerical results 
for finite plate - inertial flow with experi- 


ments (Rq/L = 3/4). 
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Figure 35. - Schematic of 
steady state normal gravity 
impingement. 
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Figure 34. - Photograph of sharp edge disks. 





czzle diameter, 0.5 cm; jet velocity, 286 cr 'sec 










Figure 36. - Continued. 




(c) Disk radius, 2.0 cm. 
Figure 36. ~ Continued. 






Figure 36. - Continued. 



Figure 36. - Continued. 




If) Disk radius, 1.5 cm. 
Figure 36. - Concluded. 
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Figure 37. - Variation of plume width with Weber number R JR„ = 4. 3%8 + 0. 1334 We. 
















<b) Laboratory equipment. 
Figure B2, - Concluded. 





